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SURVEY 

Since  writing  my  Survey  for  the  last  Progress  Report,  I have  spent  a month 
in  Japan,  attending  the  International  Conference  on  Theoretical  Physics.  This  v/as  a 
very  interesting  meeting,  though  strenuous.  Its  most  interesting  features,  from  the 
point  of  view  of  the  Group,  were  the  symposium  on  molecular  theory,  held  at  Nikko, 
in  the  mountains,  before  the  main  conference;  and  the  opportunity  to  visit  the  research 
group  in  molecular  theory  headed  by  Professor  Kotani  of  the  University  of  Tokyo. 

These  meetings  brought  together  Coulson,  Lowdin,  Mulliken,  and  a number  of  other 
workers  in  the  field,  and  gave  an  excellent  opportunity  to  talk  over  many  of  the  prob- 
lems on  which  we  are  working  with  other  workers  in  the  field,  particularly  the  Japan- 
ese. Their  programs  are  very  similar  to  ours;  and  in  some  problems,  such  as  the 
oxygen  molecule,  they  are  extremely  close.  Thus,  Ishiguro  is  carrying  through  a 
configuration  interaction  treatment  of  oxygen,  using  the  same  configurations  which 
MecVl^r  had  used,  but  with  much  better  atomic  orbitals;  his  results  so  far,  though 
very  preliminary,  look  rather  similar  to  Meckler's.  It  was  t itremely  interesting  to 
talk  over  such  questions  as  the  computation  of  molecular  integrals,  with  representa- 
tives of  most  of  the  groups  actively  working  on  them. 

Dr.  Lowdin,  who  was  in  Japan,  has  now  joined  our  Group,  for  several  months, 
and  has  already  been  giving  inter  esting  seminar  talks,  both  on  the  determination  of 
analytic  atomic  wave  functions,  about  which  he  writes  in  this  Report,  and  on  general 
methods  of  numerical  integration  of  differential  equations. 

Perhaps  the  most  interesting  new  development  since  the  last  Progress  Re- 
port has  been  the  discovery  of  the  new  method  of  handling  the  difference  equations  en- 
countered in  the  theory  of  perturbed  wave  functions  in  crystals,  which  is  reported  in 
contributions  Nos.  1 and  2 to  this  Report.  This  method,  which  was  largely  the  work 
of  Dr.  Koster,  is  much  more  elegant  and  general  than  the  method  mentioned  in  the 
preceding  Progress  Report,  and  should  open  up  the  field  of  wave  functions  and  energy 
levels  of  crystals  containing  impurity  atoms  in  a way  which  may  well  lead  to  substan- 
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me  difference  equations  will  be  met  in 


other  connections.  Thus,  the  two-electron  example  of  ferromagnetism,  about  which 
Drs.  Koster  and  Statz  and  I wrote  recently,  encounters  just  such  mathematical  prob- 
lems, and  these  new  methods  should  help  greatly  in  carrying  that  theory  further. 

Then  again,  the  exciton  problem  is  another  case  where  similar  methods  are  applicable. 
Putting  the  matter  a little  more  generally,  the  whole  question  of  electron  correlation 
is  essentially  a matter  of  the  relative  motion  of  electron  pairs  under  the  action  of 
their  Coulomb  interactions,  and  of  the  periodic  potential  in  which  they  move  in  a real 
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crystal;  and  the  separation  of  variables  in  such  a problem  into  coordinates  of  the  cen- 
ter of  mass,  and  relative  coordinates,  will  in  every  case  lead  to  a problem  of  this 
type,  so  that  we  should  be  able  to  make  real  progress  toward  a theory  of  correlation 
energy  in  periodic  potentials,  rather  than  being  limited  to  free  electron  gases  as  in 
most  previous  treatments. 

The  calculation  of  energy  bands  in  solids  is  a problem  in  which  we  are  al- 
ways very  interested,  partly  on  account  of  treating  perfect  crystals,  and  partly  be- 
cause its  results  are  necessary  in  order  to  handle  the  questions  of  impurity  atoms, 
ferromagnetism,  and  correlation  energy  about  which  we  have  just  been  speaking.  The 
augmented  plane  wave  method,  about  which  we  wrote  in  the  preceding  Progress  Re- 
port, still  seems  like  a promising  method  of  attack,  and  Dr.  Howarth  reports  on  the 
work  which  he  is  starting,  to  program  this  method  for  the  Whirlwind  digital  computer. 
As  soon  as  the  difficulties  (purely  technical)  of  this  programming  are  ironed  out,  we 
hope  to  give  the  method  a real  test,  on  copper,  and  see  whether  it  is  really  any  better 
than  the  other  existing  methods.  In  connection  with  this  method,  or  any  method  of 
energy  band  computation,  one  needs  information  about  the  atomic  wave  functions  and 
potentials,  and  for  this  reason  the  type  of  study  which  Dr.  Lowdin  is  giving  to  inter- 
polation methods  of  deriving  atomic  wave  functions  is  of  great  promise. 

In  molecular  problems,  we  continue  to  feel  that  the  most  pressing  question  is 
to  investigate  methods  of  simplifying  configuration  interactions,  so  that  we  can  handle 
more  complicated  molecules  by  methods  which  are  likely  to  be  valid.  The  study  which 
Dr.  McWeeny  is  making  of  configuration  interaction  in  simple  molecules,  proceeding 
by  the  use  of  orthogonal  atomic  crbitals,  is  likely  to  be  very  useful  in  this  respect. 

His  discovery  that  the  exchange  integrals  in  this  method  are  almost  negligibly  small 
may  prove  to  be  a real  key  to  the  simplification  of  the  calculation.  The  small  size  of 
these  integrals  was  already  known  for  the  hydrogen  molecule,  but  it  now  appears  to 
be  general  enough  so  that  perhaps  it  always  holds.  In  addition  to  McWeeny's  study, 
the  work  of  Kaplan,  and  of  Schweinler  and  Koster,  continues;  Callen's  problem  is 
practically  completed.  Meckler's  mechanization  of  the  molecular  orbital  method  con- 
tinues Pratt's  further  study  of  antiferromagnetism  demonstrates  how  really  compli- 
cated that  problem  is  In  the  preceding  Progress  Report,  I mentioned  that  something 
like  the  superexchange  which  has  so  often  been  postulated  for  this  problem  seemed  to 
form  the  basis  of  the  action  of  oxygen  in  antiferromagnetic  crystals,  and  no  reason 
has  appeared  to  change  this  conclusion.  However,  when  we  try  to  extend  Pratt's  re- 
sults to  more  complicated  cases  than  the  three-atom  example  which  he  discussed,  we 
meet  difficulties  in  the  way  of  handling  the  required  configuration  interactions  which 
are  of  the  same  sort  which  we  encounter  in  any  complicated  molecular  problem.  It  is 
clear  that  the  simple  sort  of  way  out  of  these  difficulties  which  is  usually  considered 
in  the  literature  is  by  no  means  adequate;  in  fact,  if  the  various  calculations  described 
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qualitatively  in  the  literature  are  carried  out  quantitatively,  it  seems  clear  that  they 
result  in  ferromagnetic,  rather  than  antiferromagnetic,  interactions.  We  shall  ap- 
parently have  to  go  further  in  understanding  configuration  interaction  before  we  find 
really  adequate  techniques  for  solving  the  problem  of  antiferromagnetism. 

The  various  other  subjects  under  study  by  members  of  the  Group  are  pro- 
gressing; and  with  the  addition  of  the  foreign  visitors  who  are  here  at  present,  there 
is  a lively  atmosphere  which  should  lead  to  rapid  progress.  As  mentioned  in  the  pre- 
ceding Progress  Report,  I sun  lecturing  this  term  on  energy  levels  in  perturbed  crys- 
tals, and  the  notes  of  these  lectures  will  appear  as  Technical  Report  No.  5,  which  we 
hope  will  appear  at  approximately  the  same  time  as  this  Progress  Report. 

J.  C.  Slater 
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1.  WAVE  FUNCTIONS  FOR  IMPURITY  LEVELS 


In  the  preceding  Progress  Report,  we  outlined  a general  approach  to  the 
problem  of  setting  up  wave  functions  for  a crystal  perturbed  by  impurity  atoms,  writ- 
ing the  perturbed  wave  function  as  a superposition  of  Wannier  functions  in  the  well- 
known  way-  We  proposed  that  instead  of  regarding  the  coefficients  U (R.)  of  the  Wan- 
nier  functions  of  the  n n band  on  the  atomic  site  at  vector  position  R.  as  a continuous 
function  of  a vector  position  r,  taking  on  the  values  R.  at  the  atomic  sites,  and  using 
the  approximate  differential  equation  for  the  continuous  function,  we  might  better  op- 
erate with  the  coefficients  directly,  setting  up  the  difference  equations  which  deter- 
mine them.  We  mentioned  that  in  some  important  cases  these  difference  equations 
could  be  handled  rigorously  Since  writing  that  Report,  we  have  found  a much  more 
powerful  method  of  handling  the  difference  equations,  which  is  so  useful  that  it  indi- 
cates that  this  technique  may  become  really  usable  for  practical  problems.  This 
method  will  be  outlined  in  the  present  note,  and  its  application  to  a three-dimensional 
problem  will  be  presented  in  the  next  section. 

We  start  our  procedure  by  expanding  the  perturbed  function,  not  in  terms  of 
Wannier  functions  directly,  but  in  terms  of  Bloch  functions;  though  later  we  transform 
it  into  an  expansion  in  terms  of  Wannier  functions,  which  is  physically  more  inform- 
ing. We  shall  formulate  it  first  for  the  simplest  case,  of  a single  non-overlapping 
band  with  a perturbation  only  at  the  atom  at  the  origin,  but  then  shall  describe  its 
extension  to  a general  case.  Let  us  use  periodic  boundary  conditions,  with  N atoms 
in  the  repeating  region.  We  have  Wannier  functions  a(r  - R.)  for  the  band  in  question. 

T’  in  the  corresp  jnding  Bloch  functions  are 

B^(r)  = N'  1/2  £(R.)  e*  ' ^ a(r  - R.).  (l.  l) 

These  functions  are  normalized  and  orthogonal,  and  there  are  N of  them  correspond- 
ing to  Nk  values.  The  matrix  component  of  the  perturbative  energy  between  two  of 
these  Bloch  functions  will  be 


N 


1 £(fl.  R ) el("  * ' **i  + vUl,  R.), 

J J 


(1-2) 


where  V(lt.,  R.)  is  the  matrix  component  of  the  perturbation  between  the  corresponding 


Wannier  functions 


in  our  special  case,  we  assume  that  V(lt.,  R.)  is  zero  unless  both 


R.  and  R^  are  zero,  in  which  case  it  is  VQ.  Then  the  matrix  component  (1.2),  between 
the  Bloch  functions  corresponding  to  propagation  constants  k and  k\  is  merely  V /N, 
the  same  for  all  pairs  of  propagation  constants.  It  is  this  very  simple  form  of  the 
matrix  component  which  leads  to  a very  simple  solution  for  this  problem. 

Let  us  now  set  up  our  solution  as  a series  of  Bloch  functions,  £(k  ) f^»Bj^(r ). 
Then,  remembering  that  the  periodic  potential  problem  has  a diagonal  matrix  with 
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respect  to  the  Bloch  functions,  with  matrix  components  Ej^,  the  equations  for  the  fj^-'s 
become 


EiT fiT  + I(k')(vo/N)  f-,  = Ef-.  (1.3) 

We  can  immediately  solve  these  equations,  in  the  form 


(Vq/N)  Z(iT')f-, 


(1.  4) 


Since  Z(k')  fjj*i  a constant  independent  of  k,  this  at  once  tells  us  that  f£-  is  a constant 
divided  by  E - E^*-,  which  gives  a complete  solution  of  the  problem  of  determining  the 
wave  function;  the  constant  can  be  determined  by  the  normalization  condition.  We  can 
at  once  find  an  equation  for  the  energy.  Let  (Vq/N)  Z(k')  fj^i  = A.  Then  from  (l.  4), 
fj^-  = A/(E  - E^*-).  We  substitute  this  in  the  definition  of  A above,  and  find  that 

(vo/N>  Z(k ) e4ej.  = 1 °'5) 

This  is  an  equation  for  E,  determining  the  energies  of  the  perturbed  states  as  a func- 
tion of  the  perturbative  potential  V . 

The  general  nature  of  Eq.  (l.  5)  allows  us  to  make  some  deductions  about  the 
energy  levels.  The  function  Z^  ) 1 / (E  - E^*-)  becomes  infinite  when  E equals  any  one 
of  the  Ej^-'s  Its  reciprocal,  which  according  to  (l.  5)  equals  V /N,  becomes  zero  at 
each  of  the  E^-'s,  so  that  if  Vq  is  zero,  the  solutions  are  just  the  Ej^-'s.  Between  each 
pair  of  zeros  of  this  reciprocal  function,  there  will  be  a simple  pole.  Thus  as  Vq 
goes  from  zero  to  positive  or  negative  infinity,  each  of  the  energy  levels  will  be  dis- 
placed, but  not  further  than  the  next  adjacent  E^-.  There  are  exceptional  cases,  where 
the  energy  lies  outside  the  unperturbed  energy  band.  These  are  the  discrete  levels 
produced  by  impurity  atoms.  As  tne  energy  gets  well  outside  the  unperturbed  band, 
we  may  approximate  each  term  E - Ej^  by  E - E , where  E is  the  average  energy  in 
the  band.  Thus  the  summation  will  be  approximately  N/(E  - Eq),  so  that  (1.5)  will 
become  approximately  E = E + V . A discrete  energy  level,  in  other  words,  is 

GO 

pushed  out  of  the  band,  which  becomes  linear  in  VQ  fordarge  perturbations. 

We  can  find  the  wave  functions,  as  well  as  the  energy.  We  know  that  the 
wave  function  can  be  written  as  Z(k  ) fj^Bj^r  )•  where  fj^  = A/(E  - E^-).  We  are  more 
interested  in  the  wave  function  as  an  expansion  in  Wannier  functions,  rather  than  in 
Bloch  functions.  We  can  find  this  at  once,  by  substituting  for  the  Bloch  functions  in 
terms  of  the  Wannier  functions,  from  Eq.  (l.  l).  When  we  do  this,  we  find  that  the 
coefficient  U(R^)  of  the  Wannier  function  a(r  - R.)  is  given  by 


-5- 


(WAVE  FUNCTIONS  FOR  IMPURITY  LEVELS 


U(R  ) = AN'  1/2  I(iT)  ^ (I.  6) 

This  simple  result  is  the  complete  solution  of  the  problem  of  finding  the  wave  function 
in  the  presence  of  an  impurity  producing  a perturbation  on  only  one  atom.  It  is  to  be 
understood,  of  course,  that  the  E which  appears  in  (1.6)  is  that  which  results  from  the 
solution  of  (l . 5). 

We  shall  now  state  the  results  which  hold  in  the  more  complicated  case  where 
the  perturbative  potential  is  felt  over  a number  of  atoms,  and  where  a number  of  over- 
lapping bands,  with  their  corresponding  Wannier  functions,  are  concerned;  though  we 
shall  not  give  the  analytical  treatment  here  We  find  that  the  solution  can  be  written 
as  a linear  combination  of  functions  of  the  type  we  have  just  been  considering,  arising 
from  perturbations  on  single  atoms.  One  such  function  is  to  be  centered  on  each  atom 
on  which  there  is  a perturbative  potential,  and  one  is  to  be  constructed  from  the  Wan- 
nier function  of  each  band  concerned.  The  number  of  coefficients  of  these  functions, 
then,  will  equal  the  number  of  atoms  on  which  there  is  perturbative  potential,  multi- 
plied by  the  number  of  Wannier  functions  concerned  on  each  atom.  These  coefficients 
are  determined  by  a set  of  simultaneous  homogeneous  linear  equations,  which  are  in- 
compatible unless  the  determinant  of  coefficients  vanishes.  This  furnishes  a secular 
equation  for  the  energy,  though  it  is  not  of  the  usual  type  of  secular  equation.  All  the 
various  functions  which  are  superposed  to  find  the  final  solution  are  to  be  computed 
for  the  energy  determined  from  this  secular  equation.  Consideration  of  symmetry 
properties  of  the  resulting  solution  will  ordinarily  lead  to  relations  between  the  coeffi- 
cients, reducing  greatly  the  order  of  the  secular  equation  to  be  solved;  in  some  cases 
it  can  remove  entirely  the  requirement  of  solving  a secular  equation,  and  lead  di- 
rectly to  the  final  wave  function.  This  general  method  has  not  yet  been  used  for  the 
solution  of  any  problems;  it  will  be  discussed  in  somewhat  more  detail  in  Technical 
Report  No-  5,  which  is  under  preparation. 

Reference 
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lecular Theory  Group,  M.  I.  T.  , October  15,  1953,  p 6. 

J.  C.  Slater  and  G.  F-  Koster 
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2 AN  EXAMPLE  OF  IMPURITY  LEVELS  IN  SOLIDS 


The  general  method  of  treating  impurity  levels  in  solids  mentioned  in  the  pre- 
ceding contribution  to  this  Quarterly  Progress  Report  can  be  carried  quite  far  in  the 
case  of  the  simplified  three-dimensional  impurity  problem.  Let  us  imagine  a simple 
cubic  lattice  whose  primitive  translations  are 

Kpqr  = Pal  + ^2  + ra3 

where  a^  and  a^  form  the  sides  of  the  unit  cell  and  lie  along  the  x,  y,  and  z direc- 
tions respectively.  In  this  solid,  we  could  discuss  the  effect  of  a perturbation  on  a 
single  non-overlapping  energy  band.  The  Bloch  functions  of  this  band  we  shall  denote 
by  Bj^-(r  ) and  are  eigenstates  of  a one-electron  periodic  Hamiltonian  H with  eigenvalues 


HB-(r)  = E-  B-(r)  (2.1) 

This  energy  band  is  describable  in  terms  of  a Wannier  function  a(r  - R ) which  has 
cubic  symmetry.  For  simplicity  we  shall  assume  that  these  Wannier  functions  have 
only  nearest  neighbor  interactions.  In  this  case,  the  energy  band  energies  are  given 

by 


Ek  = 2^(1)  [cos  + 


cos  k,  + cos  k. 


(2.  2) 


Here 


g{\)  = J r(T  - Rooq)  H a(r  - RlQQ)  dr 

= y*a( r - Rqoo)  K a(r"-  RQ10)  dr,  etc. 


(2-  3) 


k.,  k^,  and  k^  are  the  x,  y,  and  z components  of  the  propagation  vector  k times  the 
interatomic  distance.  We  now  apply  a perturbation  V(r)  to  this  problem  which  has  the 
property  that 


J a^r  ' K000^  H a^r  ‘ K000^  dT  _ Vo 

The  matrix  elements  of  this  perturbation  between  the  Wannier  functions  vanish  unless 
both  Wannier  functions  are  centered  about  the  atom  at  the  origin. 

From  the  preceding  report,  we  see  that  if  we  try  to  expand  our  impurity  wave 
function  in  terms  of  the  Bloch  functions  for  the  unperturbed  problem 
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(AN  EXAMPLE  OF  IMPURITY  LEVELS  IN  SOLIDS) 


Z(k)  BiT(^ 


(2.  5) 


then  for  this  case  we  have  as  the  eigenvalue  condition  that 


1 = I (Vq/n)  Me  - 2^(l)(cos  kj  + cos  + cos  k.j)J| 
kik2k3 


(2.6) 


If  the  wave  functions  of  the  perturbed  problem  are  expressed  in  terms  of  the  coeffi- 
cients of  the  Wannier  functions  at  the  various  atomic  sites  (Up^),  these  are  given  by 


U 


pqr 


= const.  Y.  exp(ikjp  + ik^q  + ik^r)/  |e  - 2<^(  1)  |cos  k^  + cos  k.,  + cos  k3J| 


k,k2k3 


(2.  7) 


If  we  allow  the  crystal  to  become  infinite,  we  can  convert  the  summations  in  (2.  6)  and 
(2.  7)  to  integrations  over  continuous  variables  in  the  first  Brillouin  zone.  In  this  case, 
the  eigenvalue  condition  becomes 


1 = V0'/w^y*  J f /e1  - (cos  kj  + cos  k^  + cos  k3jl  1 dkj  dk2  dk3  (2.  8) 

o *o  *o  ^ *■ 


The  coefficients  of  the  Wannier  functions  are  given  by 


U = const 
pqr 


it  .it  .tt  cos  kjP  cos  k^q  cos  k^r 


• Iff 


J r 

o^o  E1  - cos  k j 


+ cos  k^  + cos  k3 


dkj  dk^  dk3 


(2.9) 


Here 


E'  = E/ (l)  Vo'  = Vo/2<f(l)  (2.10) 

In  order  to  study  the  solutions  of  this  problem  all  we  need  do  is  evaluate  the 
integrals  (2.  8)  and  (2.  9).  We  shall  do  this  for  the  bound  states  in  the  case  that  Vq  and 
$ (l)  are  both  negative  (The  bottom  of  the  band  is  E 1 = 3.  ) The  integrals  (2.  8)  and 
(2.  9)  can  be  evaluated  by  doing  two  of  the  integrations  directly,  but  this  becomes  in- 
volved with  complete  ellipti  c integrals  and  the  final  integration  becomes  a difficult 
numerical  integration.  Instead  of  this  we  shall  digress  on  a different  approach  to  this 
integral. 

In  the  integral  (2.  9),  we  can  carry  out  a power  series  expansion  of  the  de- 
nominator (since  E'  > 3)  and  get  the  result  that 
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mu  I 

(dk)  cos  pk1  cos  qk^  cos  rk^  £ (coskj  + cosk,  + cosk,)n/E  n 

kj  kj  yj  n 

(2.11) 

(p  = q = r = 0 gives  the  eigenvalue  condition  integral.  ) Introducing  a parameter  t we 
see  that 


U = (1/E')  V — 
pqr  ' L n 


",n  jr  jt  pir  t(cos  kj  + cos  + cos  k^ 

— / / / cos  pk.  cos  qk,  cos  rk,  e 
dtnJoJo-lo  1 2 3 


(F>£  (l/E'"{S>(t)  Iq(,)Ir(t> 


J=0 
(2.  12) 


t = 0 


j»TT  t C O k 

I (t)  = (l/v)  I e cos  pk  dk  is  the  Bessel  function  of  pure  imaginary  argument. 

Ir"  9^ 


We  can  formally  write 

U = (tt3/E')  Tl/(  1 — )I  (t)  I (t)  I (t)]t 

pqr  > L 1 E'  dt  P q rwJt  = i 


(2-  13) 


(We  get  Eq.  (2.  12)  by  expanding  the  denominator.  ) If  we  define 

U (t)  = (tt3/E ’)  fl/(l  - — — ) I (t)  I (t)  I (t)l 
pqr  ' L ' E'  dt  P w qw  rwJ 


(2.  14) 


j 

for  all  t,  we  can  formally  multiply  both  sides  of  this  equation  by  (E‘ ) and  find  that 

dt 


(E1  - A)  U (t)  = ir3fl  (t)  1 (t)  I (t)l 
dt  pqr'  ' L P q r'  ’} 

We  can  solve  this  first  order  differential  equation  and  get  the  result  that 


(2.  15) 


- E’t  TI 
e U 


pqr 


(t)  = it3  f e E 1 I (t)  I (t)  I (t)  dt  + const. 
* * pqr 


(2.  16) 


The  constant  of  integration  can  be  shown  to  be  zero  by  studying  the  asymptotic  forms 

of  (2.  12)  and  (2.  16)  for  large  t.  Therefore  the  quantity  we  want,  U (0)  = U , is 

M J ' pqr'  pqr’ 

given  by 


fe'E,t  Ijt)  IJt)  I„(t)  dt 


P q r' 


(2.  17) 


The  eigenvalue  condition  is  given  by 
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(AN  EXAMPLE  OF  IMPURITY  LEVELS  IN  SOLIDS) 


(2. 


The  wave  function  (2.  17)  is  not  normalized  but  the  normalization  can  easily  be  found 


l(k)  flfr  = £(k)(E  - E^f  2 = - ± I(iT)(E  - E^)-  1 =J ^ te'  E,t  [l0(t)j3  dt 

r/*00  -E't  r ■) - 1/ 2 

Therefore  (2.  17)  can  be  normalized  by  multiplying  by-J  / t e Iq(i.)  dt> 

After  this  digression  on  the  evaluation  of  this  integral,  we  can  give  the  re- 
sults on  the  values  of  the  wave  function  and  energy.  The  integral  in  (2.  18)  was  car- 
ried out  for  a number  of  values  of  E'  and  the  corresponding  value  of  the  perturbation 
V ' was  found.  These  results  are  given  in  Fig.  2-1.  We  notice  that  a finite  value  of 
the  perturbation  is  necessary  before  the  bound  state  pulls  out  of  the  bottom  of  the  en- 
ergy band.  For  V ' slightly  greater  than  this  minimum  perturbation  the  energy  be- 
haves quadratically  as  a function  of  V '.  For  large  perturbations  the  energy  E'  ap- 
proaches V ' and  becomes  essentially  linear  as  a function  of  V '. 

The  coefficients  of  the  Wannier  functions  have  also  been  computed  for  the 
energies  in  Fig.  2-1.  These  coefficients  have  been  computed  (unnormalized)  in  the 
(100)  directic..  and  are  given  in  Table  2-1.  The  values  in  the  expression  (2.  17) 

uPoo  - fo  ''"''[yfyi d* 


Table  2-1 


Up00  - f e E t Ip(t)  [Iq^)]2  dt  as  a function  of  energy  and  p 


X 

3.  0 

3.  1 

3.  2 

3.  5 

4.  2 

5.  0 

0 

. 4990 

.4311 

. 4001 

.3410 

. 2646 

. 2137 

i 

. 1658 

. 1123 

. 0935 

. 0645 

. 0371 

. 0238 

2 

. 0778 

. 037b 

. 0270 

. 0 i 4 i 

. 00564 

a a 

. UU£  i O 

3 

. 0501 

.0154 

. 00929 

. 00350 

. 000924 

. 000343 

4 

. 0411 

. 00717 

. 00365 

. 000964 

. 000162 

. 0000439 

It  has  already  been  pointed  out 


(1) 


that  if  the  difference  equations  for  the  Wan 


nier  function  coefficients  were  replaced  by  an  approximate  differential  equation  we 
might  expect  the  coefficients  of  the  Wannier  functions  to  fall  off  as  e ^r/r  (r  f 0) 


(AN  EXAMPLE  OF  IMPURITY  LEVELS  IN  SOLIDS) 


Fig.  2-1 

Energy  (E')  vs  perturbation  potential  ( V *) 
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Fig.  2-2 

Coefficients  of  the  Wannier  function  in  the  (100)  direction 
(p(Up0Q/U  10q))  as  a function  of  p for  various  values  of  E'. 


(AN  EXAMPLE  OF  IMPURITY  LEVELS  IN  SOLIDS) 


where  y = 0 when  the  bound  state  first  comes  out  of  the  energy  band  (E1  = 3).  This 
behavior  is  also  predicted  when  we  assume  that  E*  is  near  to  the  bottom  of  the  band 
and  the  principle  contribution  to  the  integral  (2.  9)  comes  from  F^-  near  the  bottom  of 
the  band  when  the  denominator  of  the  integral  is  small.  In  this  case, 

/i(pk.  + qk?  + rk  ) 

-£ - - L—  (2.  20) 

AE!  + + k22  + k32) 

where  AE'  = E'  - 3 and  the  denominator  is  obtained  by  expansion  of  cos  k^,  etc. 
Carrying  out  the  integral  (2.  20)  gives  the  result  that 

U 

pqr 

which  agrees  with  the  result  of  the  differential  equation  approach. 

In  Fig.  2-2  we  plot  p(up00/U)00)  on  semi -log  graph  paper.  The  log  of  this 
function  should  be  a linear  function  of  p if  the  behavior  predicted  by  the  differential 
equation  approach  or  the  approximation  to  the  integral  (2.  9)  is  valid.  We  have  not 
plotted  the  points  Uqqq  since  neither  approximate  approach  is  valid  for  this  point. 

We  notice  that  for  E'  = 3 the  graph  is  very  rearly  linear  and  horizontal.  This  is  pre- 
dicted by  both  approximations.  For  E1  > 3 the  linearity  becomes  quite  good  for  points 
removed  from  the  impurity  atom.  (i.  e.  , p = 2,  3,  4,  lie  nearly  on  a straight  line. ) 

We  can  also  compare  the  actual  slopes  of  the  curves  in  Fig.  2-2  with  those  predicted 
from  approximation  (2.  21).  Eq.  (2.  21)  predicts  a slope 


For  E'  = 3.  1 this  gives  a value  . 193,  for  E'  = 3.  2 the  value  is  . 275.  The  correspond- 
ing slopes  calculated  from  the  actual  values  of  U^qq  are  . 209  for  E'  = 3.  1 and  . 285 
for  E'  = 3.  2.  We  can  see  therefore  that  the  approximation  to  the  slope  gives  a value 
to  within  about  5%  of  the  true  value. 

Work  on  this  problem  will  continue  to  calculate  the  coefficients  of  the  Wan- 
jnier  functions  in  other  directions  than  the  ( 1 0 0 ) direction,  in  order  to  check  the 
spherical  symmetry  of  the  wave  function. 

Reference 

1.  J.  C.  Slater  and  G.  F.  Koster,  Quarterly  Progress  Report,  Solid-State  and  Mo- 
lecular Theory  Group,  M.I.  T.  , October  15,  1953,  p.  7. 

G.  F.  Koster 
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3.  ENERGY  BANDS  IN  COPPER 


The  use  of  the  recently  proposed  augmented  plane  wave  method^  for  calculat- 
ing energy  bands  in  crystals  depends,  for  its  practical  use,  upon  the  convergence  of 
the  secular  equation  which  arises  as  an  attempt  to  ensure  continuity  of  slope  of  the 
wave  function  at  the  boundary  of  the  matching  sphere.  It  is  hence  of  immediate  im- 
portance to  apply  the  method  to  a particular  metal  to  study  the  convergence  and  value 
of  the  method.  The  case  of  copper  has  been  selected  as  being  of  particular  interest 
owing  to  the  overlapping  of  the  4s  conduction  band  with  the  narrow  3d  band.  Calcula- 
tions of  the  3d  band  structure  of  copper  have  been  made  previously  by  the  tight  binding 
method,  ^ and  of  the  3d  and  4s  band  by  the  cellular  method  ^ There  is  close  re- 
semblance in  the  results  for  the  d band  in  these  two  cases,  and  it  is  hoped  that  the 
present  method  will  agree  substantially  with  these  calculations,  and  in  addition  may  be 
used  to  investigate  the  higher  excited  electronic  states. 

The  initial  step  in  the  augmented  plane  wave  method  is  to  calculate  the  energy 
E of  the  augmented  wave  in  terms  of  the  energy  Eq  of  the  plane  wave  outside  the  match- 
ing sphere  of  radius  R The  dependence  of  E upon  Eq  is  given  by 


F(E)  = «(E  - Eq) 


4w  X (21  + 1)  j/U/EoR)  K^ln  Uj(E,  r) 


r = R 


(3.  1) 


U is  the  volume  of  the  unit  cell  outside  the  matching  sphere,  which,  for  a face-centered 
cubic  structure  such  as  copper  is 


Q = 


a 

1 


3 


a being  the  lattice  constant.  jj(x)  is  the  spherical  Bessel  function,  and  Uj(E,  r)  is  the 
radial  wave  function.  Both  the  numbers  of  parameters  involved  and  the  fact  that 
Uj(E,  r)  is  obtained  by  numerical  integration  have  led  to  an  attempt  being  made  to  com- 
pute (3.  1)  by  use  of  the  Whirlwind  high  speed  digital  computer. 

The  problem  falls  into  three  portions:  the  computation  of  Pj(E,  r),  defined  by 

i'Uj(E,  r)  = P|(E,  r); 


the  evaluation  from  inis  of 


Vj(E,  r)  ~ rZ  In  u^E,  r) 

and  lastly,  the  computation  of  (3.  1)  as  a function  E for  various  values  of  Eq,  a,  R. 

The  program  devised  will  compute  (3.  l)  and  determine  its  roots  for  a given  number  of 
lattice  parameters  a,  for  a given  number  of  matching  radii  R for  each  a,  thus  enabling 
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us  to  see  the  effect  on  the  energy  levels  of  choosing  different  matching  radii,  and  for  a 
given  number  of  equally  spaced  energies  Eq. 

The  function  P^E,  r)  is  a solution  of  the  equation 

= f(p,  r)  = P fv  - E + + il]  (3.  2) 

dr^  - r^  J 


where  V is  the  spherically  symmetrical  potential  in  which  the  electron  is  supposed  to 
move  inside  the  matching  sphere.  In  order  to  forward  integrate  (3.  2),  use  is  made  of 
the  Gauss  Jackson  formula 


P 


n 


12  n 


+ h 


(3.3) 


where  6 2 fn  represents  the  second  sum  of  f and  h is  the  interval  of  r.  The  first 
correction  term  which  has  been  omitted  in  (3.  3)  is  - h2/240  6^  fn>  which  can  be  made 
negligible  by  choice  of  sufficiently  small  intervals  h.  In  order  to  use  (3.  3)  to  advance 
the  integration  from  rn  ; to  r^,  the  second  sum  6 2 f is  computed  from  previous 
knowledge,  and  is  given  as 


6' 2 f = 26*  2 f 

n n - 1 


f 

n - 


1 


In  general,  the  resulting  algebraic  equation  (3.  3)  must  be  solved  for  P using  (3.  2), 
by  iteration.  In  the  case  of  the  linear  equation,  P^  can  be  determined  directly,  the 
method  reducing  to  the  Noumerov  method. 

To  preserve  generality,  the  general  solution  of  (3.  2)  involving  iterations  has 
been  set  up  on  Whirlwind.  It  is  planned  to  increase  the  speed  by  taking  into  account 
the  linearity  of  (3.  2),  to  avoid  the  iterations. 

The  subroutine  to  advance  the  integration  by  one  step  is  incorporated  in  a 
master  routine  which  organizes  the  change  of  interval  at  given  points,  and  the  storage 
of  the  wave  function.  In  order  to  carry  out  the  rest  of  the  computation,  we  require  a 
table  of  the  wave  function  for  values  of  r greater  than  a prescribed  value  Rq.  To  use 
this  program,  the  data  to  be  specified,  which  is  dependent  upon  the  par  ticular  metal 
under  consideration,  consists  of  the  potential  V(r)  tabulated  at  intervals  of  r,  in,  in 
the  range 


r. 

i 


^ r £ r{ 
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where  r =0,  and  i <:  k,  all  of  which  parameters  must  also  be  specified.  In  the  case 

of  copper,  we  are  using  k = 2,  and  h,  = . 005,  h,  = . 04,  r = . 04,  r,  = 4.  0.  In 

(0)  1 ( i ) c 1 c (c) 

addition,  two  starting  values  Pj  ',  Pj  ',  together  with  the  initial  value  of  r,  r'  ’ are 

required  to  start  the  integration.  The  next  section  of  the  program  computes  v^(E,  R) 
for  a given  E and  f,  and  for  various  R.  A table  of  the  derivatives  of  P is  computed 
from  the  table  of  P by  a forward  difference  formula.  For  a given  R,  P and  its  deriva- 
tive are  interpolated  from  these  tables  by  use  of  the  Everett  modified  second  difference 
formula,  using  three  entries  on  either  side  of  the  argument;  the  truncation  error  is 
about  0.  0005  6^P.  The  function  v^(E,  R)  may  thus  be  computed  and  stored  for  all  R.; 
the  preceding  computation  is  repeated  for  all  the  i values  required,  and  lastly,  the 
parameter  E is  increased  by  equal  amounts  between  two  specified  values.  It  is 
planned  to  use  a sufficiently  small  range  of  E for  each  set  of  integrations  carried  out 
to  enable  the  same  starting  solutions  P^\  P^1^  to  be  used  for  all  E.  The  particular 
radii  R and  the  number  of  parameters  i,  E,  are  the  only  data  required  to  carry  through 
this  part  of  the  program. 

Lastly,  the  use  of  v^(E,  R)  to  compute  (3.  1)  in  the  form  Farj,  (E)  is  a straight- 
forward procedure.  The  only  operation  requiring  comment  is  the  calculation  of  the 
spherical  Bessel  functions  j^(  v^E^K).  These  are  computed  from  the  slowly  varying 
function  Aj+  defined  by 


For  high  i,  Aj  + , j ^ can  be  calculated  from  a polynomial  expansion,  and  use  is  then 
made  of  the  backward  recurrence  formula 


2 

AJ- l/2(x)  = AJ+i/2(x)  ‘ (fj-  f)(zi  +~)  Ai+  3/2W 

This  proceeds  with  negligible  loss  of  accuracy  for  all  f when  x < 10,  which  is  sufficient 
for  the  present  purpose.  The  convergence  of  the  sum  in  (3.  1)  is  governed  by  the  con- 
vergence of  jj(  v^qR);  for  values  of  Eq  corresponding  to  wave  vectors  spanning  the 

first  Brillouin  zone,  terms  for  which  J > 6 are  neelipible.  while  for  lartrer  E . it  is 

w “ “ o 

estimated  that  terms  with  i > 12  may  be  neglected.  The  first  two  portions  of  this  pro- 
gram have  been  tested  and  made  error  free.  The  last  section,  together  with  methods 
for  determining  the  roots  of  (3.  1)  is  at  present  being  tested,  preparatory  to  using  the 
complete  program  to  investigate  the  energy  band  structure  of  copper. 
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4.  THE  AUGMENTED  PLANE  WAVE  METHOD 


I have  begun  work  to  apply  the  augmented  plane  wave  method,  described  in 
the  previous  Progress  Report,  to  sodium.  An  easy  generalization  of  that  method  has 
been  made  which  may  or  may  not  prove  to  be  useful.  The  generalization  is  this: 
ra.her  than  choose  our  wave  function  to  be  continuous  across  the  matching  spheres, 
we  choose  the  derivative  of  the  wave  function  to  be  continuous  across  them.  This 
condition  along  with  the  minimization  of  the  expectation  of  the  energy  relative  to  the 
energy  inside  the  spheres  yields  wave  functions  with  properties  analogous  to  the  wave 
functions  generated  by  the  original  procedure.  All  functions  with  different  energy 
but  with  the  same  k value  are  orthogonal  and  have  no  component  of  the  energy  between 
them. 

The  E vs  k1*  curve  for  these  functions  has  asymptotes  for  values  of  E for 

— ♦ 

which  u^'(E;  R.)  vanishes.  The  branches  of  these  curves  touch  the  asymptotes  for  k 
values  which  cause  the  derivative  of  a spherical  Bessel  function  j.'(kR.),  to  vanish. 
Thus  if  we  superimpose  the  E vs  k curves  for  the  two  cases  (continuous  slope  method 
and  continuous  wave  function  method)  we  find  that  between  every  two  old  asymptotes 
we  find  a new  one,  and  between  every  two  old  points  of  tangency  to  the  asymptote  we 
find  a new  one.  Furthermore,  the  slope  of  each  branch  is  everywhere  decreasing. 
These  comments  are  illustrated  in  Fig.  4-1. 


Fig.  4-1 

A typical  branch  of  E vs  in  the  old  and  new  schemes 

These  new  functions  may  be  useful  in  speeding  up  the  convergence  of  the 
secular  equation  in  that  their  spectrum  combined  with  the  old  spectrum  yields  a den- 
ser spectrum  than  either  alone.  This  affords  us  more  functions  for  our  secular  equa- 
tion which  have  energies  close  to  the  energy  we  wish  to  calculate.  It  may,  however, 
not  pay  to  make  this  combination  but  to  use  one  scheme  for  one  range  of  energies  and 
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the  other  scheme  for  another  range.  These  ideas  will  be  investigated  in  the  applica- 
tion of  the  method  to  sodium. 

Because  the  asymptotes  are  suggestive  of  "gaps"  it  was  thought  that  the 
scheme  could  be  modified  to  place  these  "gaps"  at  more  natural  energies  than  those 
which  cause  Uj'(E;R.)  or  u^(E;R.)  to  vanish.  The  modification  was  to  make  the  com- 
bination, a(wave  function)  + p(slope  of  wave  function),  (where  a and  p are  constants 
to  be  determined)  continuous  across  the  boundary  spheres.  The  expectation  of  energy 
of  such  a wave  function  is  then  minimized  with  respect  to  the  energy  interior  to  the 
spheres  and  the  parameter  p/a.  This  yields  two  equations.  We  can  then  solve  for 
P/a  as  a function  of  energy  and  then  solve  for  energy  as  a function  of  k^.  The  "gaps" 
will  then  occur  for  values  of  energy  which  satisfy 


ui'(E;Ri) 

uJ(E;Ri) 


£ = y(E) 


The  "gaps"  will  no  longer  occur  at  the  same  energy  for  all  k but  actually  change  with 
k.  Another  attractive  feature  is  that,  since  we  have  minimized  the  expectation  of 
energy  relative  to  two  parameters,  our  approximation  to  the  actual  E vs  k.  curves 
should  be  better  than  in  the  two  special  cases  in  which  we  take  either  a = 0 or  p = 0 
to  start  with.  However,  the  mathematical  difficulties  in  solving  the  energy  equation 
for  this  scheme  are  so  formidable  that  I do  not  consider  its  application  worth  while  at 
the  present  time. 


M.  M.  Saffren 
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5.  ENERGY  BANDS  IN  CHROMIUM 


The  calculation  of  matrix  elements  of  the  Hamiltonian  between  two  Bloch 
waves  (3d  or  ion-core)  is  in  progress.  The  calculation  of  matrix  elements  of  this 
type  has  been  completed  for  the  short-range  portion  of  the  crystal  potential.  (The 
short-range  portion  of  the  crystal  potential  is  a spacial  sum  over  the  crystal  lattice 
of  non-overlapping  atomic  potentials;  the  remainder  of  the  crystal  potential,  composed 
of  a spacial  sum  of  strongly  overlapping  atomic  potentials,  is  called  the  long-range 
portion.  ) Work  is  continuing  on  the  calculation  of  matrix  elements  of  the  kinetic  en- 
ergy and  the  long-range  portion  of  the  crystal  potential. 

R.  H.  Parmenter 
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6.  STUDIES  OF  ATOMIC  SELF-CONSISTENT  FIELDS 


The  study  of  the  properties  of  the  atomic  self-consistent  fields,  which  was 
actually  started  some  time  ago  at  Uppsala,  has  been  continued  during  my  stay  at  the 
Solid-State  and  Molecular  Theory  Group.  Of  the  main  problems  to  be  solved,  we  have 
treated  the  question  of  the  analytic  representation  of  the  self-consistent  field  (SCF) 
wave  functions,  the  interpolation  (or  extrapolation)  of  these  wave  functions  to  atoms 
or  ions  for  which  they  have  not  yet  been  computed,  and  the  problem  of  the  connection 
between  the  SCF -functions  and  the  simple  hydrogen-like  ones. 

The  question  of  the  analytic  representation  has  been  solved  in  two  steps.  The 
numerically  given  SCF-functions  were  first  expressed  in  a form,  originally  proposed 
by  Slater,  ^ as  hydrogen-like  functions  having  each  exponential  replaced  by  a sum  of 
one,  two,  three,  or  more  exponentials.  The  exponents  and  coefficients  of  these  ex- 
ponentials have  been  determined  by  a method  of  "quotient  analysis"  and  successive 
approximations,  and  the  details  and  the  results  of  this  first  attack  have  already  been 
published.  ^ However,  the  method  is  tedious  and  rather  time-consuming,  since  the 
time  required  for  an  almost  perfect  fit  (0.  001)  is  usually  of  the  order  30-40  hours  by 
using  an  ordinary  desk  machine. 

The  method  for  finding  the  analytic  fit  has  now  been  modified.  In  the  new 

approach,  the  SCF -function  is  expressed  as  the  product  of  a polynomial,  having  the 

same  nodes  r r ".  . . , as  the  wave  function  itself,  and  a sum  of  exponentials.  The 
o o 

essential  problem  is  therefore  to  develop  the  function 

g(r)  = fni(SCF  r)/rJ+  1(rQ'  - r)(rQ"  - r)  . . . , (6.  l) 

taken  in  equidistant  points,  into  a sum  of  geometrical  series,  each  one  corresponding 
to  an  exponential  term.  In  principle,  such  a problem  can  be  exactly  solved,  but  due  to 
the  vanishing  of  valid  figures,  this  exact  solution  is  without  practical  value.  Instead 
we  have  used  a method  of  successive  approximation  starting  from  the  tail  of  the  func- 
tion, where,  in  addition  to  the  previously  mentioned  quotient  analysis,  we  have  also 
used  the  theory  for  developing  a series  of  equidistant  entries  into  two  geometrical  pro- 
gressions. The  technique  is  so  simple  and  straightforward  that  it  now  takes  about  4- 
5 hours  to  get  a perfect  analytic  fit  even  for  a rather  extensive  SCF-function  having 
several  nodes. 

The  analytic  SCF-functions  have  been  utilized  for  a first  attack  on  the  inter- 

( 1 ) 

polation  problem.  We  have  retained  the  rule  given  by  Slater  ' that  the  exponents 
should  be  linear  functions  of  the  atomic  numbers,  but  we  have  replaced  his  rule  for 
the  coefficients  by  a simpler  law  based  on  the  analogy  with  the  hydrogen-like  case. 

Test  interpolations^  using  these  rules  have  shown  an  almost  perfect  accuracy  (0.  001  - 

0.  002). 

We  have,  however,  also  attacked  the  interpolation  problem  in  a more  direct 
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way  by  using  only  the  numerical  tables  of  the  SCF-functions.  We  have  used  the  basic 
principle  that  any  quantity,  which  for  a hydrogen-like  function  would  equal  the  atomic 
number  Z,  for  a SCF-function  corresponds  to  an  "effective  nuclear  charge"  Ze^,  and 
practically,  we  have  then  found  that  this  Zg^  is  an  almost  linear  function  of  the  atomic 
number  well  suited  for  interpolation  purposes.  By  considering  the  corresponding  Z^^, 
we  have  in  this  way  treated  the  interpolation  of  the  positions  and  magnitudes  of  the 
maxima  and  minima,  the  derivatives  in  origin,  the  eigenvalues,  etc. 

An  important  point  is  that  the  positions  of  the  eventual  nodes  may  be  inter- 
polated with  almost  perfect  accuracy  in  this  way.  The  nodes  usually  cause  complica- 
tions in  the  problem  of  interpolating  the  SCF-functions  as  a whole,  and  we  have  there- 
fore found  it  convenient  to  divide  them  away  and  tc  .onsider  simply  the  function  g(r) 
given  by  (6.  1).  The  logarithm  of  g(r)  is  for  a hydrogen-like  function  a linear  function 
of  rZ,  and,  for  a SCF-function,  the  corresponding  quantity  is  therefore  a linear  func- 
tion of  rZ£^(r)  where,  for  a given  r,  the  second  factor  is  an  almost  linear  function  of 
the  atomic  number.  The  interpolation  problem  may  therefore  simply  be  solved  by 
making  a linear  interpolation  in  the  logarithm  of  g(r)  for  each  r-value  under  considera- 
tion with  an  eventual  renormalization  of  the  interpolated  wave  function  as  a whole. 

Test  interpolations  have  shown  that  this  simple  rule  gives  an  accuracy  of  0.  001  for  the 
inner  wave  functions  and  an  accuracy  of  about  0.  005  or  better  for  the  outermost  one. 

In  this  way  it  is  also  possible  to  compute  an  effective  charge  Zg^(r)  which  is 
a continuously  varying  function  of  the  radius.  In  addition  to  this  function,  we  have 
studied  the  quotient 

«*>(r)  = Zeff(r)/Zeff(0).  (6.2) 

which  is  monotonously  decreasing  from  the  value  w(C)  = 1.  Practically  we  have  found  that 
this  function  is  very  slowly  varying  with  the  atomic  number,  and  it  may  therefore  be 
used  for  improving  the  accuracy  of  the  interpolation  of  the  outermost  function.  Another 
characteristic  feature  is  that,  for  an  isoelectronic  series,  the  difference  between  the 
two  u-functions  for  the  fields  with  and  without  exchange  seems  to  be  approximately  the 
same  function  of  r for  all  elements  of  the  series.  This  property  may  be  used  for  es- 
timating wave  functions  with  exchange  from  those  without  exchange,  when  only  one 
field  with  exchange  is  known  in  the  series.  For  the  moment,  several  fields  with  ex- 
change have  been  determined  up  to  Cu+  (Z  = 29)  whereas  various  fields  without  ex- 
change have  been  computed  up  to  Hg  (Z  = 80). 

As  applications  of  the  interpolation  methods  here  described  we  have  estimated 

a 3 

"initial"  wave  functions  for  two  Ne-like  configurations,  namely  Ne  and  Al'  , and 
carried  through  the  integrations  of  the  Hartree-Fock  equations  for  fields  with  exchange. 
The  corrections  to  the  "final"  wave  functions  were  extremely  small. 


-22- 


(STUDIES  OF  ATOMIC  SELF-CONSISTENT  FIELDS) 


Finally,  we  have  investigated  the  closer  connection  between  the  SCF-functions 
and  the  hydrogen-like  functions  by  using  a method  with  a varying  scale  factor.  A func- 
tion F(r)  is  here  said  to  be  obtained  from  another  function  f(r)  by  means  of  a constant 
scale  factor  X,  if 


F(r)  = X1/2  f(Xr),  (6.  3) 

which  transformation  leaves  the  normalization  integral  invariant.  An  obvious  generali- 
zation of  this  relation  is  to  replace  the  linear  function  Xr  by  a more  general  function 
ili  = 4*(r ) and  to  put 

F(r)  = [^]i/2  f^),  (6-4) 

which  transformation  also  leaves  the  normalization  integral  invariant.  At  the  same 
time  we  obtain 


^ F2(|)  d£  = d£, 

# o •'o 


(6.  5) 


which  relation  shows  that  we  simply  can  get  a table  of  the  function  i|i(r)  by  inverse  in- 
terpolation from  tables  of  the  total  spherical  charges  corresponding  to  the  two  radial 
distributions  under  consideration.  The  work  is  essentially  simplified  by  preparing 
the  table  for  the  right-hand  member  of  (6.  5)  with  a particularly  small  interval,  but  this 
work  may  fortunately  be  done  just  once  for  all  if  we  use  the  same  set  of  standard  func- 
tions f(r),  for  instance  the  hydrogen  functions. 

According  to  this  method  we  are  able  to  express  all  SCF-functions  with  and 
without  exchange  as  hydrogen-like  functions  having  an  effective  charge  given  by 


zeff(r)  = +(r)/r. 


(6.6) 


All  Is-,  2s-,  and  2p-functions  available  have  been  investigated  in  this  way,  and  we 
have  found  that  even  the  ^-functions  are  very  convenient  for  interpolation  purposes. 


TV, 


v*  method  is  simple  and  pu.tc  straightforw*ard,  but  it  shows  a complication  fur  me 


neighborhood  of  the  nodes  which  appear  to  be  slightly  "out  of  scale".  The  method 
gives  also  a possibility  of  deriving  the  mutual  scaling  of  two  SCF-functions  on  each 
other. 


The  effective  charge  functions  derived  in  this  way  may  be  useful  for  estimat- 
ing molecular  integrals  for  SCF-functions  from  the  standard  tables  for  hydrogen-like 
functions  given  by  Kotani,  Kopineck,  and  others.  Let  us  consider  a diatomic  molecule 
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of  two  identical  atoms  with  the  internuclear  distance  equal  to  a.  The  molecular  tables 
are  then  usually  given  in  terms  of  an  effective  nuclear  charge  Z,  and  we  have  now  at- 
tempted to  get  the  SCF-value  simply  by  replacing  this  Z by 

Z = Zeff(rQ  + a/2),  (6.  7) 

where  rQ  is  a constant  determined  by  the  condition  that  the  approximation  should  be 
valid  for  a = 0.  For  the  2pn-functions  of  carbon,  which  are  particularly  important  in 
the  theory  of  conjugated  organic  compounds,  we  have  carried  out  a comparison  between 
these  estimations  and  the  actual  integral  values  obtained  by  numerical  integration:  the 
accuracy  of  the  estimations  appeared  to  be  surprisingly  high,  mostly  around  95  - 100%. 
However,  for  a-functions,  the  method  must  certainly  be  further  developed  before  a 
similar  accuracy  can  be  obtained. 

All  the  numerical  computations  are  now  essentially  concluded  ana  a series  of 
articles  is  being  prepared  for  publication. 
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7.  A REVISED  VALENCE  BOND  METHOD 


Work  has  proceeded  along  the  lines  indicated  in  the  previous  Report  (to  be 
referred  to  as  I).  Numerical  calculations  on  the  four -electron  problem  with  the  sym- 
metry of  the  square  (in  fact,  the  v-electron  system  of  the  hypothetical  molecule  cyclo- 
butadiene) have  been  completed  and  the  results  are  sufficiently  interesting  to  warrant 
inclusion  in  the  present  Report.  More  interesting  will  be  the  results  from  the  six- 
electron  problem  (the  hexagonal  benzene  molecule):  calculations  are  already  well- 
advanced  and  it  therefore  seems  worth  while  reporting  first  on  the  general  VB  approach 
which  is  common  to  both  problems. 

In  the  four-electron  case  the  following  VB  structures  (i.  e.  , spin-paired  or 
valence  bond  structures  defined  over  strictly  orthogonal  orbitals,  AO's)  are  employed 
(see  I for  further  details) 

Singlet  States 

a<-ba  b a *-  b a b 

I t 

d -*•  c d c d -*  c d -*  c 

+ - - + 

2 non-polar  8 adjacent-charge  polar 
Triplet  States 


b 

X.  - 


+ 

a b 


+ 

a 


4 nonadjacent -charge  polar 


d -»  c d -*■  c 
+ “ - + 
6 doubly-polar 


a — b a — b a — b a b 

\4  ... 

d-c  d — c d - c d/Xc 

+ ~ ~ + 

3 non-polar  8 adjacent-charge  polar  4 nonadjacent -charge  polar 

The  structures  not  shown  are  obtained  simply  by  rotation  of  the  given  pictures  and  in 
each  case  the  basis  of  structures  is  complete  subject  to  the  limitation  imposed  by  us- 
ing a basic  set  of  2pTr  orbitals  alone.  In  the  triplet  case  there  are,  at  first  sight,  four 
non-polar  structures;  but  only  three  are  linearly  independent,  for  a "vector  rule"  for 
the  arrows  shows  that  the  fourth  is  (-1)  times  the  sum  of  the  other  three.  The  num- 
ber of  linearly  independent  spin  eigenfunctions  of  total  spin  S and  given  Sz  which  can 
arise  from  a given  orbital  configuration,  and  therefore  the  number  of  structures  show- 
ing a given  dispostion  of  charges,  follows  of  course  from  the  Wigner  formula  as 
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where  N is  the  number  of  singly  occupied  orbitals;  this  provides  a simple  and  useful 
check  on  the  number  of  structures  set  up  by  pictorial  methods. 

In  the  six-electron  problem  a little  consideration  shows  that  there  are  175 
independent  singlet  structures  and  189  triplets;  and,  although  the  hexagonal  symmetry 
permits  a very  substantial  reduction  of  tne  problem,  use  of  the  full  number  of  con- 
figurations would  incur  both  the  evaluation  of  many  hundreds  of  matrix  elements  and 
the  solution  of  secular  equations  of  rather  high  order  (e.  g.  , 22  singlet  combinations 
of  Ajg  symmetry  can  be  formed).  There  are  clearly  two  reasons  for  wishing  to  se- 
lect, from  the  formidable  array,  just  those  structures  of  greatest  potential  value  in 
the  energy  calculation;  in  the  first  place,  extension  of  these  methods  to  even  slightly 
more  difficult  problems  will  be  quite  impossible  unless  a useful  method  of  selection 
can  be  found;  and  secondly,  our  physical  insight  into  the  calculation  is  progressively 
weakened  as  more  and  more  configurations  are  admitted.  The  methods  now  being  de- 
veloped provide  a very  natural  and  plausible  way  of  making  the  necessary  selection: 
for  (see  I)  non-polar  and  singly-poiar  structures  are  expected  to  predominate  in  the 
description  of  a bonded  system  --  pictorially,  they  are  sufficient  to  describe  "charge 
hops”  across  all  possible  bonds  and  the  inclusion  of  multiply-polar  structures  (which 
are  in  any  case  energetically  less  well  favored)  could  only  fulfil  somewhat  more  freely 
this  fundamental  bonding  requirement.  First  then  we  should  look  for  the  best  approxi- 
mation which  can  be  built  from  non-polar  and  singly-poiar  structures;  this  would  re- 
duce the  problem  to  one  involving  only  65  singlet  and  99  triplet  structures.  We  notice 
however  that  all  bonds  between  nearest  neighbors  can  be  accounted  for,  pictorially,  in 
terms  of  non-polar  and  ad jacent -charge  polar  structures;  in  a preliminary  calculation 
then  we  might  tentatively  admit  this  additional  restriction  --  working  in  terms  of  29 
singlet  and  45  triplet  structures.  This  last  step  in  the  selection  becomes  fairly  con- 
vincing when  it  is  realized  that  admission  of  nonad jacent -charge  structures  is  tanta- 
mount to  a recognition  of  considerable  interactions  between  atoms  which  are  known  to 
be  "non-bonded".  At  the  same  time  the  conventional  VB  structures,  which  arise  by 
simple  generalization  of  the  Heitler-London  type  function  and  which  presumably  do 
in  some  way  describe  the  bond  systems  to  which  they  are  supposed  to  correspond,  may 
be  expanded  in  terms  of  VB  structures  and  the  importance  of  a structure  is  then  found 
to  fall  off  as  the  separation  cf  the  charges  increases. 

The  types  of  structure  we  shall  first  consider  are  shown  below  (together  with 
the  number  of  structures  of  each  type  --  formed  by  rotating  the  pictures) 
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Singlet  Structures 
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(3) 


non-polar 


a-±-b 


c f 


\ ✓ 

e d 


(6) 


(45  in  all) 


a--b 


t f 


a b 


a — b 


a — b 


c f c f c f- 

\ ✓ \ ✓ \ ✓ 
e d e d e d 

(6)  (6)  +(6) 

adjacent-charge,  singly  polar 


■ c f ■ 


e d 

(6)+ 


a-±-b 

e <-  d 

(6) 


It  will  be  notice  that  crossed-link  structures  are  linearly  dependent  on  those  given  above, 
as  also  are  types  such  as 


which  is  equivalent  to 


Rules  for  picking  a linearly  independent  set  are  given  elsewhere.  ^ 

It  is  now  worth  discussing  some  aspects  of  the  actual  calculations.  The  vari- 
ous structures  may,  in  both  problems,  be  formed  into  combinations  providing  bases  for 
irreducible  representations  of  the  point  group  of  the  molecule:  this  part  of  the  calcula- 
tion is  quite  straightforward  and  will  not  be  described  here  --  but  it  will  be  noted  that 
the  structures  possess  particularly  simple  transformation  properties,  merely  being 
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permuted  under  the  group  operations  (a  property  not  shared  by,  for  instance,  MO  con- 
figurations). The  matrix  eigenvalue  equation  then  breaks  into  "block  form",  each 
block  involving  structure-combinations  of  like  symmetry  and  having  elements  which 
may  be  reduced,  using  the  group  properties,  into  expressions  which  contain  a rela- 
tively small  number  of  distinct  elements  of  H between  pairs  of  individual  structures, 
and  the  latter  may  be  evaluated  without  difficulty  by  the  rule  given  in  the  last  Report. 
The  application  of  this  rule  (which  was  quite  general)  to  VB  structures  will  now  be  il- 
lustrated by  reducing  two  typical  matrix  elements  completely  to  one-  and  two-electron 
integrals.  Take  the  pairs  of  structures  below;  they  may  be  described  as  showing,  re- 
spectively, a shift  of  negative  charge  and  a shift  of  positive  charge,  and  each  is  a trip- 
let structure. 


a — b 


f 

\ 

e 


c 

✓ 

d 


+ 


a — b 


f c 


A B 


C D 


(1)  (AIHIB) 

The  cycle  diagram  is  (ignoring  at  first  the  arrow  directions) 


and,  by  direct  application  of  the  rule,  gives  a matrix  element 

22|Qd(ce)  - + Kdb^ce^  ' (Kdf(ce^  ‘ 1 Kdf(re^] 

where  the  modified  Coulomb  and  exchange  integrals  are,  more  fully  (see  also  I)  Q^(ce)  = 
(clhie)  + (aclglae)  + (bclglbe)  + ...  = Pce  + (aa,  ce)  + (bb,  ce)  + ...  and,  for  instance, 
Kda(ce)  = (ca|g|ae)  = (ca,  ae).  Here  Pcg  is  the  "resonance"  integral  of  an  MO  treat- 
ment and  vvc  use  a common  notation  in  which  (ab,  cd)  is  the  mutual  potential  energy  of 
charge  distributions  ab  and  cd.  In  the  present  case  obvious  symmetries  exist,  it  is 
convenient  to  employ  certain  standard  Coulomb  integrals  --  such  as  Q (ab),  Q (ac),  and 
the  simplest  form  of  the  result  is 

[(Q a(aC)  + (aa>  ac)  " (bb,  ac))  - ~(ac,ae)  - j(ab,ad)j. 
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Since  the  original  pictures  are  "matched"  by  an  even  number  of  arrow  reversals, 
there  is  no  need  to  affix  a factor  - 1 and  this  is  the  final  result.  In  practice  the  re- 
duction described  is  immediate. 

(2)  (CIHID) 

The  cycle  diagram  is 


but  in  this  case  only  a one -orbital  difference  really  occurs,  for  the  columns  containing 
orbitals  d and  c in  the  determinants  of  D may  be  interchanged,  a factor  - 1 being 
simultaneously  added;  and  it  is  evident  that  this  preliminary  matching  is  precisely 
equivalent  to  making  the  interchange  in  the  cycle  diagram  and  breaking  the  cycle  to 
give 


(the  (-)  indicating  that  the  matrix  element  is  to  be  prefixed  by  a factor  - 1).  The  rule 
is  now  directly  applicable,  giving  the  result 

(CIHID)  = -21[Qce(dd.ec)  - ±(2Ked(ec))  + Kef(ec)  - Kea(ec)  - Keb(ec)J 

which  again  reduces  on  inspection  to 

-2  [(Q  a(ac)  + (ac,  bb)  - (ac,  aa))  - (ab,  be)  - (ac,  ae)J 

In  working  with  valence  bond  structures  a cycle  change  of  this  type  is  appropriate 
whenever  a shift  of  positive  charge  is  involved,  but  in  no  other  case. 

Although  all  the  matrix  elements  involved  may  be  written  down  with  relatively 
little  labor,  by  the  above  methods,  it  must  not  be  forgotten  that  the  integrals  which 
occur  are  defined  over  the  orthogonal  AO's  (though  we  shall  now  use  the  bar  notation 
only  when  actually  making  comparisons  between  VB  and  VB  quantities).  In  fact,  all 
the  overlap  and  many-center  integrals  which  the  conventional  theory  omits  are  im- 
plicitly included  (but  in  a more  manageable  way)  in  the  simple  types  of  AO  integrals 
which  now  appear.  Thus  the  distinct  AO  integrals  of  the  type  (pq,  rs),  which  may  be 
set  out  as  a matrix  (A)  with  elements  labeled  by  two  double  suffixes,  A^  rs(P  < q,  r < s), 
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transform  into  integrals  (pq,  rs)  when  the  AO's  are  replaced  by  AO's.  And  if  the  or- 
bital  transformation  be  written  (p)  = (p)(\)  then  (A)  = (\  ) (A)(\  )g  where  (\  )g  is  the 

symmetrized  Kronecker  square  of  (\).  This  transformation  (which  in  the  six-electron 
case  involves  a product  of  three  21  X 21  matrices  can  obviously  be  reduced  when  mo- 
lecular symmetries  exist  --  but  not  in  general.  It  may  therefore  ultimately  be  worth 
programming  this  transformation,  which  would  also  describe  the  calculation  of  inte- 
grals over  MO's,  for  digital  computation.  Fortunately,  however,  when  the  new  or- 
bitals are  strongly  localized,  an  integral  (pq,  rs)  is  quite  negligible  unless  the  atom 
pairs  pq,  rs  are  very  near  (excluding  p = q,  r = s)  and  then  depends  only  on  their  im- 
mediate environment:  thus,  in  practice,  the  AO  products  such  as  pq  may  be  expanded 
directly,  only  a few  neighboring  orbitals  occurring  with  significant  weight,  the  im- 
portant terms  in  the  integrand  of  (pq,  rs)  then  being  picked  out  by  inspection. 

Certain  features  of  the  transformation  to  AO's  have  been  noted  in  the  previous 
Report.  Broadly  speaking,  integrals  of  the  form  (pp,  qq)  do  not  differ  greatly  from 
their  conventional  counterparts  (pp,  qq)  --  but  all  other  integrals  (of  types  (pq,  rr)  and 
(pq,  rs))  are  almost  completely  "transformed  away".  This  reduction  in  numerical 
values  of  all  the  exchange  and  hybrid  integrals  is  so  striking  that  it  is  worth  giving 
figures  (see  Table  7- 1 ):  the  AO  values  which  follow  are  taken  from  Craig^  in  the  first 
case  and  from  Parr,  Craig  and  Ross^  in  the  second  (all  values  are  in  e.  v.  ).  Most 
of  the  integrals  not  given  are  of  the  order  0.  01  or  less;  and  there  is  clearly  a strong 
case  for  neglecting,  in  AO  calculations,  most  of  the  troublesome  many-center  terms 
which  cannot  be  discarded  in  other  formulations.  T'  ere  are,  of  course,  reasons  for 
expecting  an  effect  of  this  kind  --  for  a charge  distribution  pq  contains  zero  net  charge 
and  interactions  between  differently  located  distributions  should  therefore  be  of  "sec- 
ond order"  --  but  the  results  more  than  fulfill  our  expectations:  indeed,  it  seems 
likely  that  the  complete  neglect  of  all  non-Coulombic  (or,  more  loosely,  "exchange" 

--  including  hybrid)  integrals  may  have  little  effect  upon  our  ultimate  general  conclu- 
sions. * We  shall  therefore  go  on  to  examine  the  energy  calculation,  verifying  previous 
conclusions^  about  the  form  of  the  secular  problem  (see  I),  and  investigating  the  ef- 
fect of  the  many  small  exchange  integrals. 

Below  is  given  the  H -matrix  (which  is  symmetric),  referred  to  the  basis  of 


*It  will  be  noticed  that  the  assumptions  made  by  Pariser  and  Parr^/  in  a recent  serni- 
empirical  MO  theory,  which  at  first  sight  seem  unreasonably  drastic,  may  now  be 
simply  interpreted:  indeed,  if  the  AO's  are  used  in  constructing  a MO  theory,  our  ob- 
servations on  the  numerical  values  of  the  integrals  encountered  suggest  immediately 
the  very  assumptions  made  by  Pariser  and  Parr.  The  success  of  our  scheme  of  neglect- 
ing all  exchange  integrals  would  thus  simultaneously  justify,  in  a non-empirical  way, 
the  general  sufficiency  of  the  framework  developed  by  those  writers  (essentially  an 
MO  theory),  at  the  same  time  providing  a more  satisfactory  theoretical  interpretation 
of  the  basic  parameters  appearing  in  their  formulation  as  integrals  over  AO's. 
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5 Aj  singlet  functions.  The  types  of  structure  from  which  these  functions  are  con- 
structed are  indicated  by  the  labels  on  the  rows  and  columns:  n and  p signify  non- 
polar and  polar  structures,  affixes  1 or  2 meaning  singly  or  doubly  polar  and  suffixes 
a or  na  meaning  adjacent  or  non-adjacent  charges. 


Table  7-2 


n 


na 


n 

P 

P 

P 

P 


(Q-  0,244)  2(Pab+  0.058)  2 + 0.  284) 


1 

a 

1 

na| 

2 

a 

2 


(Q+  9.  063) 


ac 


- Ji( 0.  280) 

(Q  + 10.  895) 


*a 

2(0.  039) 
2(pab+  0.058) 

2v^((3ac  + °-284) 

(Q+  21.  842) 


2^(0.  026) 
2V/2(P  i t 0.  058)1 

3D  j 

0 

2 72(0.  122) 

(Q  + 15.  166) 


The  Q appearing  along  the  diagonal  is  appropriate  to  a non-polar  configuration,  the 
much  higher  energy  Q-ierms  for  the  polar  structures  being  expressed  as  Q plus  the 
increase  in  Coulomb  repulsion  energy  accompanying  the  charge  transfer;  thus,  for  an 
adjacent  charge  structure  the  Q-term  is  Q (bb)  = Q + (aa,  aa)  - (aa,  bb)  =•  Q + 

17.  537  - 8.  790.  The  difference  between  this  value  and  the  diagonal  element  shown 
arises  from  the  small  exchange  terms  which  come  along  with  Q (bb).  The  only  other 
large  elements  are  those  containing  the  p's  (Craig  uses  p , = - 2.857,  p = - 0.  355) 
these  arising  between  those  structure -combinations  which  contain  members  differing 
by  a shift  of  one  charge.  It  is  convenient  to  refer  all  energies  to  Q as  zero,  writing 
x = E - Q,  and  with  this  notation  the  eigenvalue  equation  is  (H)(c)  = x(S)(c),  where, 
in  this  instance  (with  normalized  orthogonal  combinations),  (S)  is  the  unit  matrix.  Nu- 
merical values  of  the  elements  of  (H)  are  given  below  with  and  without  the  exchange 
terms:  and  it  is  obvious  that  the  general  pattern  of  the  secular  problem  is  not  depend- 
ent on  their  inclusion. 


Table  7-3 


0 - 5.  72 

- 1.  02 

0 

0 

- 0.244  - 5. 

598 

- 0.  201 

0.  078 

0.  074 

8.  75 

0 

- 5.  72 

-8.  09 

9. 

063 

- 0. 396 

- 5.  598 

- 7.  917 

11.  01 

1.  02 
21. 98 

0 

0 

15.  24 

10.  895 

0.  201 
21. 842 

0 

0.  345 
15.  166 

Approximate  Accurate 
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Calculations  have  been  carried  out  for  the  lowest  singlet  states  of  symmetry 
Aj  and  and  for  the  lowest  triplet  of  symmetry  A^g.  The  approximate  eigenvec- 
tors and  energies  (relative  to  Q)  are  shown  below. 


Table  7-4 


n 

Pa 

p 

na 

Pa 

Pa 

Accurate 

0.  806 

0.  537 

0.  051 

0.  Ill 

0.  219 

Approximate 

0.  802 

0.  538 

0.  051 

0.  120 

0.  224 

n 

Pa 

Pa 

Accurate 

0.  821 

0.  563 

0.  099 

x = - 6. 426 

Approximate 

0.  822 

0.  557 

0.  118 

X 

II 

1 

O' 

*0 

Ni 

n 

Pa 

Prta 

Accurate 

0.  776 

0.  576 

0.  274 

x = - 5.  915 

Approximate 

0.  773 

0.  574 

0.  273 

x = - 6. 00 

x = - 3.  871 
x = -4.  06 


The  most  striking  conclusion  to  be  drawn  from  these  results  is  that  the  exchange  in- 
tegrals do,  indeed,  have  no  very  great  effect  on  our  predictions.  If  this  is  generally 
valid  (and  calculations  on  benzene  so  far  provide  no  evidence  to  the  contrary)  an  im- 
mediate simplification  of  the  whole  problem  can  be  effected:  for,  ignoring  from  the 
start  the  many  small  exchange  terms,  matrix  elements  vanish  for  all  pairs  of  struc- 
tures which  correspond  to  more  than  one  charge  shift.  And,  for  normalized  structures, 
the  non-vanishing  elements  are  given  by  the  extremely  simple  rules 

Table  7-5 


*r*2 

Non-polars 

Similarly  charged 
polars,  a+  b- 

Non-polar  with 
polar  a+  b" 

Polar  a*  b- 
with  a+  c " 

Polar  a"  b+ 
with  a"  c + 

(+1ihi+2) 

2n'nOQ 

2n‘n°  Qa(bb) 

^n'"°eab 

2.  2n  ~ n°  p 

'DC 

- 2.  2n  "no  p 

rbc 

33- 


(A  REVISED  VALENCE  BOND  METHOD) 


where  n is  the  number  of  cycles  in  the  diagram  for  4^  and  (and  is  easily  visualized 
as  the  number  of  "islands"  arising  from  superposition  of  their  link  patterns)  and  nQ 
is  the  total  number  of  links  in  any  structure  (=  N / Z);  these  results  are,  as  usual,  to 
be  prefixed  by  - 1 if  spin-matching  involves  an  odd  number  of  arrow  reversals.  Thus, 
in  this  approximation,  all  the  matrix  elements  can  be  written  down  immediately,  with 
no  formal  analysis  whatever. 

There  is  only  one  disconcerting  feature  in  the  results  so  far:  the  doubly-polar 
structures  --  which  lie  between  20  and  30  e.  v.  above  the  ground  state  --  are  clearly 
of  considerable  importance,  and  even  if  they  are  dealt  with  by  the  approximate  method 
(i.  e.  , exchange  integrals  neglected)  the  dimensions  of  the  problem  become,  more 
generally,  depressingly  large.  First,  however,  we  shall  ask  how  good  an  approxima- 
tion can  be  constructed  from  non-polar  and  adjacent-charge  singly-polar  structures 
alone;  and  we  shall  see  how  successful  a function  is  provided  by  the  conventional  VB 
structures  --  which  can  be  expanded  in  terms  of  VB  structures^  to  give  an  approxi- 
mate eigenvector  in  the  VB  basis.  These  results  are  summarized  below. 


Table  7-6 


One  con- 

Accurate 

figuration 

MO 

Best  VB  combination 

(0. 905 

0.  424 

0 

0 

0)  x = - 2.  870 

- 3.  871 

- 0.  77 

Expanded  VB  function 

(0. 968 

0.  25? 

0 

0 

0)  x = - 2.  380 

Best  VB  combination 

(0.  837 

0.  547 

0) 

X 

= - 6.  134 

- 6. 426 

-4.  20 

Expanded  VB  function 

(0. 912 

0.  410 

0) 

X 

= - 5.  621 

Best  VB  combination 

(0.  862 

0.  508 

0) 

X 

= - 4. 773 

- 5. 915 

-4.  25 

Expanded  VB  function 

(0.  932 

0.  363 

0) 

X 

= - 4. 312 

The  energy  estimates  given  in  the  last  two  columns  are  included  for  comparison. 
These  results  prompt  two  concluding  remarks:  (1)  The  VB  structures  of  conventional 
theory  (i.  e-  , defined  over  the  overlapping  AO's)  may  be  useful  in  constructing  quite 
good  approximate  wave  functions  by  direct  expansion  into  the  more  convenient  VB 
structures.  In  the  present  case  there  is  only  one  VB  combination  of  each  symmetry, 
so  no  variational  calculation  is  required:  more  generally  the  permissible  VB  com- 
binations define  a subspace  in  the  space  of  VB  structures  and  best  weights  may  be 
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found  by  a relatively  trivial  eigenvalue  calculation.  In  any  event,  the  approximate  VB 
function  may  then  be  refined  by  an  eigenvalue  calculation  in  the  full  basis,  when  this 
is  practicable.  (2)  Even  the  crude  VB  approximation  may  be  considerably  more  ac- 
curate than  the  one-configuration  MO  approximation,  particularly  for  the  strongly 
bound  states.  Bearing  in  mind  the  example  of  the  hydrogen  molecule,  this  is  not  al- 
together unexpected;  it  appears  that  the  gross  overweighting  of  the  multiply  polar 
structures,  in  MO  calculations,  is  more  damaging  to  the  wave  function  than  their  com- 
plete neglect  (which  is  a feature  of  the  present  method  in  its  most  abbreviated  form). 
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8.  LIMITED  CONFIGURATION  INTERACTION  TREATMENT  OF  THE  NH ^ MOLECULE 

The  numerical  evaluation  of  all  necessary  integrals  for  the  observed  inter  - 
nuclear  separations  (OIS)  and  for  the  planar  model  is  still  underway.  In  the  latter 
case  it  is  assumed  that  as  the  nitrogen  atom  oscillates  from  one  side  of  the  hydrogen 
plane  to  the  other  the  N-H  distance  remains  constant.  Experimental  results  indicate 
that  this  assumption  is  quite  near  the  truth.  Computationally  this  assumption  leads  to 
desirable  duplication  of  many  radial  integrals  with  the  OIS  case. 

H.  Kaplan 
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9.  CONFIGURATION  INTERACTION  APPLIED  TO  THE  HYDROGEN  MOLECULE 


Calculations  on  this  problem  are  being  completed.  The  work  is  now  being 
written  up  as  a thesis. 


E.  C alien 


10.  THE  WATER  MOLECULE 


Since  the  last  Quarterly  Progress  Report,  we  have  obtained  and  checked  the 
expansion  for  the  charge  density  of  the  hydrogen  Is  wave  function  (orthogonalized  to 
the  oxygen  wave  functions)  about  the  oxygen  nucleus.  The  integrals  which  had  already 
been  computed  were  extensively  checked  upon  discovery  of  small  errors  in  some  of 
the  basic  functions.  These  errors  either  proved  to  have  negligible  effect  on  the  final 
result  or  have  now  been  corrected.  Work  on  the  remaining  integrals  is  continuing  as 
fast  as  our  allotted  time  on  the  I.  B.  M.  machines  permits. 

G.  F.  Koster  and  H.  C.  Schweinler 
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11.  MECHANIZATION  OF  MOLECULAR  ORBITAL  DETERMINATION 


The  machine  is  almost  ready  for  production  runs  on  Lhe  self-consistent  de- 
termination of  linear  combinations  of  orbitals.  There  remains  only  the  problem  of 
convergence  which  will  be  investigated  empirically.  The  way  the  program  is  now  set 
up,  the  machine  begins  with  a first  guess  to  the  density  matrix,  forms  the  one-electron 
effective  Hamiltonian  matrix,  diagonalizes  it,  constructs  a new  density  matrix  from 
the  eigenvectors,  computes  the  total  energy  of  the  system.  The  new  density  matrix 
is  then  fed  back  in,  unmodified,  to  begin  a new  cycle  if  the  new  and  old  total  energies 
do  not  agree  to  sufficient  accuracy.  In  the  case  of  the  Hartree-Fock  differential  equa- 
tions straight  feeding  back  of  the  newly  computed  potential  often  leads  to  a divergent 
calculation.  There  is  a good  possibility  that  our  situation  is  more  stable.  Differen- 
tial equations  correspond  to  infinite  matrical  equations,  infinite  and  continuous;  we 
deal  with  discrete,  finite  matrices.  Explicitly,  all  the  density  matrices  of  the  itera- 
tion scheme,  those  before  and  after,  differ  only  by  unitary  transformations;  the  same 
Hamilton-Cayley  equation  is  satisfied  by  all.  There  is  therefore  a certain  bounded- 
ness to  our  region  of  computation.  With  the  high-speed  machine  we  need  not  worry 
too  much  about  slow  convergence  but  only  about  out-and-out  divergence.  Suspecting 
that  divergence  is  impossible  we  will  nevertheless  try  to  force  it.  If  there  is  no  di- 
vergence we  will  learn  about  the  rate  of  convergence. 

We  will  take  Lowdin's^  exponential  fit  to  the  self-consistent  Helium  Is  func- 
tion, concoct  some  Is -like  function  orthogonal  to  it,  and  attempt  to  mix  the  two.  The 
initial  density  matrix  should  be  diagonal  --to  the  accuracy  of  Hartree  and  Lowdin. 
Insertion  of  outlandishly  wrong  first  guesses  will  yield  an  empirical  study  of  conver- 
gence. 


Reference 

1.  P. -O.  Lowdin,  Phys.  Rev.  90,  120  (j.953). 
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12.  ANTIFERROMAGNE'nSM 


In  the  previous  Progress  Report  the  details  of  the  solution  to  a three -center 
spin  coupling  problem  were  given.  An  attempt  has  been  made  to  generalize  the  re- 
sults of  this  problem  to  a larger  system  of  alternating  magnetic  and  non-magnetic 
atoms  in  a chain  or  ring.  If  such  a system  is  investigated  by  the  methods  of  configura- 
tion interaction,  the  number  of  excited  configurations  becomes  very  large  even  for  the 
simplest  model  and  this  would  lead  to  unmanageable  secular  equations.  Therefore,  it 
was  hoped  that  the  complications  inherent  in  the  rigorous  treatment  might  be  avoided 
by  using  the  results  of  the  three -center  problem  and  treating  each  magnetic  center  in 
the  chain  in  the  same  way  and  each  non-magnetic  center  in  the  same  way. 

The  following  procedure  was  used  with  these  ideas  in  mind.  The  superex- 
change configurations  of  Kramers  are  the  basis  of  the  spin  coupling  in  the  three-center 
problem.  Let  this  system  be  taken  as  Aj — B — A^  where  the  A centers  represent 
paramagnetic  ions  and  have  one  electron  associated  with  each  center  in  the  ground  con- 
figuration, and  the  B center  represents  a non-magnetic  ion  having  two  electrons  as- 
sociated with  it  in  the  same  orbital  in  the  ground  state.  Then  the  superexchange  con- 
figurations are  those  in  which  one  and  finally  both  electrons  are  transferred  from  the 
non-magnetic  ion  to  the  magnetic  ions.  The  superexchange  treatment  might  be  ap- 
proximated by  a single  configuration  if  a compound  orbital 

B = b + C(Aj  + Az) 

is  set  up  to  describe  the  electrons  on  the  non-magnetic  center,  this  orbital  being  kept 
doubly  filled.  The  orbitals  b,  Aj,  and  A^  are  orthonormal  and  associated  with  the 
various  centers  as  indicated.  The  constant  C is  to  be  determined  by  minimizing  the 
energy.  The  singlet  and  triplet  states  of  the  three -center  system  are  given  as 

= const.  |a,ci  Ba  Bp  A.,P  t A ^ p Ba  Bp  A (12.  l) 

If  one  expands  out  the  B orbitals,  a certain  amount  of  cancellation  comes  about  due  to 
the  Pauli  principle.  The  net  result  is  that  the  B functions  reduce  to  polarized  orbitals 
of  the  type  b + CA  j and  b + CA^  directed  towards  the  Aj  and  A^  centers  respectively. 
The  singlet  state  given  above  is  essentially  the  result  of  forming  a covalent  bond  be- 
tween Aj  and  the  polarized  function  b + CA^,  and  a covalent  bond  between  A^  and  the 
polarized  orbital  directed  toward  that  center.  This  type  of  approach  could  easily  be 
extended  to  a chain  or  ring  and  would  lead  to  just  one  configuration,  that  in  which  each 
A orbital  is  singly  occupied  and  each  B function  doubly  filled. 

Although  this  type  of  approach  is  appealing  it  has  a serious  fault  which  is  of 
interest.  If  the  singlet  and  triplet  states  given  in  (12.  1)  are  expanded,  we  obtain 
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4/(S  = 0)  = const.  |4>0  + + 2C2<j>2}  (12.2a) 

4*(S  = 1)  = const.  + CfjJ  (l2.  2b) 

where  in  (12.  2a)  4>0t  ♦j#  and  <|>2  represent  the  singlets  corresponding  to  the  initial, 
singly  excited,  and  doubly  excited  configurations  respectively.  In  12.  2b)  4>Q  and  <f>j 
represent  the  triplet  states  for  the  initial  and  singly  excited  configurations. 

In  the  configuration  interaction  calculation  using  orthogonal  orbitals  and  han- 
dled in  the  usual  way,  there  would  be  two  independent  constants  in  the  final  singlet  and 
one  independent  constant  in  the  triplet.  This  approximate  treatment,  which  has  only 
one  independent  constant  as  shown  in  (12.  2a)  and  (12.  2b),  can  handle  the  triplet  case 
exactly  but  not  the  singlet  case.  If  the  constant  C turned  out  to  be  small  and  if  the 
singlet  energy  were  not  very  sensitive  to  the  amount  of  the  doubly  excited  configuration 
admixed,  within  reasonable  limits,  then  the  approximate  method  might  be  expected  to 
work.  Unfortunately,  however,  the  energy  difference  between  the  singlet  and  triplet 
states  depends  critically  on  the  weight  given  the  doubly  excited  configuration.  This  can 
be  shown  graphically  as  follows.  Let  the  triplet  state  as  found  by  configuration  inie^  - 
action  be 


<MS=1)  = const.  |<|>o  + di4»i|  (12-3) 

and  let  the  singlet  state  be 

+(S  = 0)  = const.  j$0  + c1«)>1  + c2<|>2j-  (12.4) 

Now  we  can  plot  the  singlet  energy  as  a function  of  |c,|2  meaning  that  for  a given  value 

2 ^ 
of  |c2|  we  find  the  best  singlet  energy  or  best  value  of  c^  by  the  variation  method. 

This  is  shown  in  Fig.  12-1  where  the  triplet  energy,  which  is  independent  of  the  weight 
given  the  doubly  excited  configuration,  is  also  plotted.  The  fact  that  the  singlet  must 
lie  above  the  triplet  for  |c2|2  = 0 was  discussed  in  the  last  Progress  Report.  It  is  as- 
sumed that  we  are  trying  to  describe  a system  in  which  the  singlet  does  lie  below  the 
triplet  for  some  value  of  |c2|2. 

The  importance  of  highly  excited  configurations,  when  a system  is  described 
in  terms  of  orthogonal  orbitals,  was  first  shown  by  Slater  in  his  analysis  of  the  hydro- 
gen molecule.  ^ The  same  sort  of  situation  is  encountered  here  in  that  the  type  of 
spin  coupling  predicted  is  directly  connected  to  the  role  played  by  the  highest  excited 
state  of  this  system  as  shown  in  Fig.  1.  The  singlet  obtained  from  the  approximate 
approach  in  (12.  2a)  is  qf  the  form  that  the  weight  given  the  doubly  excited  configuration 
is  directly  dependent  on  the  weight  given  the  singly  excited  state.  Therefore,  the  best 
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value  of  C obtained  from  the  variational  method  may  weight 
the  doubly  excited  configuration  in  such  a manner  that  the 
|c,|*  best  singlet  will  lie  above  the  triplet.  The  calculation  was 
carried  out  for  the  problem  described  in  the  last  Progress 
Report  and  the  singlet  was  found  to  be  of  higher  energy  than 
singlet  energy  tripiet  contrary  io  the  actual  state  of  affairs.  In  view  of 
Vtriplet  energy  the  generally  unfavorable  nature  of  the  above  approximate 
method  it  seems  of  no  value  to  carry  it  further. 

Kramers^  divised  a method  of  handling  the  con- 
Fig.  12-1  figuration  interaction  problem  in  which  the  excited  configura- 

tions of  the  system  are  treated  as  a perturbation  on  the  ground  configuration.  The  un- 

N 

perturbed  ground  states  consist  of  the  2*  single  determinants  which  differ  in  the  spin 
assignments  to  the  N singly  occupied  orbitals  of  the  ground  configuration.  The  eigen- 
states of  the  system  are  found  by  diagonalizing  an  effective  Hamiltonian  but  only  with 
respect  to  the  unperturbed  ground  states.  The  secular  problem  can  be  written 


2N 

Y(U  - E6  ) C = 0,  (12.  5) 

^ . mn  mn'  n 
n = 1 


and  where  U is  the  effective  Hamiltonian 
mn 
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mn 


H + 
mn 
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H H 
ma  an 

E - 


v y Hrna  Hap  Hpn 

-HaaKE  -Hpp)  ■" 


(12.  6) 


and  a,  p,  etc.  refer  to  the  excited  states  of  the  system.  The  spin  coupling  between  the 
non-overlapping  magnetic  atoms  first  appears  in  the  third  order  term. 

As  pointed  out  above  it  appears  to  be  a rather  general  property  of  orthogonal 
orbitals  that  when  they  are  used  to  describe  a system  one  must  take  into  consideration 
highly  excited  states  in  order  that  the  results  have  any  meaning  at  all.  The  third  order 
term  in  the  effective  Hamiltonian 


I 

a / P 


H H . H„ 

JTlfl  £•  K Kn 

(E  - Hal)(E  -fl^) 


(12.  7) 


can  only  connect  the  ground  configuration  with  a doubly  excited  configuration  at  most. 
Therefore,  it  seems  unlikely  that  one  can  obtain  reliable  results  if  the  higher  order 
terms  in  the  effective  Hamiltonian  are  neglected. 

This  method  is  easily  applied  to  the  three-center  system  Aj  — B — A^.  The 
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two  unperturbed  ground  states  are 


tj  = |AjQ  ba  bp  A^PI.  t^  = |AjP  ba  bp  A^al 


(12.8) 


and  U . is 

X ■»  X 

1 c. 

= . 2dZ(J  + g) 4dhg (12.  9) 

12  (E  - EX)Z  (E  - E1)(E  - E2) 

where 

d = (AJb)  + (AjblAjAj)  + (bblbAj)  + (bA2|AjA2) 

J = (bAjlAjb)  > 0 (12.  10) 

g = (bA  J A2b) 

h = (blAj)  + 2(A,b|A1A2)  + (bA2|A2A2) 

the  orbitals  Aj,  b,  and  A2  being  orthonormal.  E^  and  E2  are  the  energies  of  the 
states  where  one  and  two  electrons  respectively  have  beeen  transferred  to  a magnetic 
atom  from  a non-magnetic  center.  The  first  term  in  Utjt^  < 0.  then  the  triplet  lies 
below  the  singlet.  Considering  only  the  first  term  we  see  that  it  is  negr  ive  definite 
since  J is  always  greater  than  g.  This  is  in  accord  with  the  fact  that  in  a treatment 
of  the  three -center  system  which  excludes  the  doubly  excited  configuration  the  singlet 
cannot  lie  below  the  triplet  regardless  of  what  orthogonal  orbitals  are  used. 

The  second  term  in  (12.  9)  must  be  large  enough  to  reverse  the  sign  of  Ut^ 
if  the  third  order  theory  is  to  be  adequate  to  describe  a three-center  problem  where 
the  singlet  does  lie  below  the  triplet.  In  order  for  the  second  term  to  make  this  change 
the  following  inequality  must  be  satisfied 

lh|  > <12  111 

Since  g is  often  quite  small  the  inequality  may  not  hold  and  one  would  be  forced  to 

carry  higher  terms  in  the  effective  Hamiltonian  even  for  this  three-center  problem. 

(4' 

The  perturbation  treatment  of  superexchange  has  been  applied  by  Shimizu'  ' 
to  a 3N  electron  system  but  restricted  to  singly  excited  configurations  only.  When  one 
uses  the  fact  that  exchange  integrals  are  never  negative  when  expressed  in  terms  of 
orthogonal  orbitals  (upon  which  this  whole  approach  is  based),  and  if  it  is  assumed 
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that  intra-atomic  exchange  integrals  are  in  general  greater  than  inter -atomic  exchange 
integrals,  then  Shimizu's  resluts  show  a ferromagnetic  spin  coupling.  This  result  al- 
so seems  to  show  the  need  for  including  more  highly  excited  states  in  the  theory.  Once 
these  states  are  admitted  the  spin  coupling  will  not  only  be  of  nearest  neighbor  type  but 
will  include  next  and  even  more  distant  neighbor  interactions.  Also  there  will  be  terms 
connecting  the  unperturbed  states  which  differ  by  more  than  the  interchange  of  a single 
pair  of  spins.  If  the  perturbation  method  is  not  rapidly  convergent,  these  new  terms 
may  well  play  an  important  role. 
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13.  NUCLEAR  ELECTRIC  QUADRUPOLE  INTERACTION  IN  THE  KC1  MOLECULE 


Further  investigation  into  the  suitability  of  the  one -electron  function, 


4-  = 4*o(l  + Crv  cos  0) 


(13.  1) 


for  describing  polarization  effects  has  been  made.  As  indicated  in  a previous  Progress 
Report,  ^ Slater  and  Kirkwood^  have  used  4< 


4-  = e r £l  + C(r  + — ) cos  ©j 


for  the  ground  state  of  the  hydrogen  atom  in  a uniform  field  and  (13.  1)  was  shown  to  be 
equivalent  to  this  in  that  they  both  gave  the  correct  value  for  the  polarizability.  Hasse^" 
also  made  very  accurate  calculations  of  the  polarizability  with  the  function 


4»  = 4> 


o |l  + (C^r  + C2r2)  cos  0J 


for  the  lowest  states  of  H and  He.  The  problem  of  the  H atom  in  a uniform  field  has 

(4) 

been  solved  exactly  by  Epstein  in  parabolic  coordinates.  Writing  his  first  order 
functions  in  spherical  coordinates  yields  a dependence 


4>  = J^l  + CR(r)  cos  oj 


for  the  excited  states  as  well  as  the  ground  state.  (R(r)  is  a polynomial  in  r for  which 

(5) 

the  highest  power  is  2 for  states  through  2p.  ) Buckingham'  ' has  obtained  good  results 
for  several  closed  shell  atoms  using  the  form  (13.  1)  with  v = i for  all  one-electron 
functions. 

For  our  F calculation  mentioned  in  the  previous  Progress  Report  we  are  thus 
led  to  the  following  ten  orthonormal  one-electron  functions: 


4*1  - N^e  r(l  + A 1 r 1 cos  0)a 

4*2  = Nje  8'75r(l+Ajr  1 cos  0)  (3 


__  : - 2.  50r 

4>3  = N2lre 


„ - 7 4«n.  . V? 

C^e  *)(1  + A^r  “ cos  9j  a 


4>4  = N2(re'2,  50r  - C1e'7'48r)(l  + A2/2  cos  9)  p 
4-5  = N3(re-2'  35r  sin  0 ei4>)(l  + cos  9)  a 

^6  = ^ 3(re  2‘  35r  sin  6 eit^))(l  + A 3r  3 cos  0)  p 


(NUCLEAR  ELECTRIC  QUADRUPOLE  INTERACTION  IN  THE  KC1  MOLECULE) 


4>7  = N3(re'  3^r  sin  6 e ^)(l  + A3r  3 cos  6)  a 

4*8  = N3(re’  35r  sin  ft  e~  ^)(l  + A^r  3 cos  8)  p 

= N^(re  3^r  - C^re  ^°r  - C^e  ^r)  CQS  Q ^ + 4 cos  Q 

4*10  = N4(re  2‘  35r  - C^re  2 50r  - C3e~  75r)  cos  0 (l  + A4r  4 cos  0)  p 

where  the  N's  represent  normalization  factors.  For  the  4*  's  we  have  used  the  functions 

(6)  ^ 

given  by  Morse,  Young,  and  Haurwitz'  ' which  have  slightly  better  values  for  the  con- 
stants appearing  in  the  exponents  than  those  given  by  Slater's  shielding  constants.  ^ 

The  above  function?  have  been  formed  into  a single  determinant  and  an  energy  and  po- 
larizability calculation  is  being  made  using  the  following  Hamiltonian: 


« = -IV  - I 


2z 

r- 

l 


I i rr 


+ <f  I r. cos  0i 


i>j 


The  A's  and  v's  are  treated  as  variation  parameters  and  orthogonality  is  maintained 
with  C j,  C^,  and  C.,.  The  integrals  for  the  angular  parts  are  obtained  using  Condor, 
and  Shortley's  C^ffm^.f  m^').  ^ The  integrals  involving  functions  of  r are  being  cal- 
culated by  methods  applied  in  a previous  Progress  Report^  and  by  use  of  the  , Lninj 

(91  * 

functions  developed  by  Schrieffer.  ' ' 
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14.  SCATTERING  OF  SLOW  NEUTRONS  BY  0? 

The  calculation  mentioned  in  the  last  Report  is  still  in  progress.  We  present 
here  an  outline  of  some  of  the  results  so  far  obtained. 

We  are  concerned  here  only  with  a neutron's  spin-spin  magnetic  interaction 
with  the  molecule,  so  that  effectively  the  neutron  interacts  only  with  the  electrons  of 
unpaired  spin.  There  are  two  of  these,  since  the  ground  electronic  level  of  O^  is  a 
triplet.  The  charge  distribution  of  the  electrons  of  unpaired  spin  has  the  same  sym- 
metry as  the  total  electronic  charge,  and  may  be  regarded  as  cylindrically  symmetri- 
cal and  fixed  relative  to  the  internuclear  axis,  since  for  slow  neutron  scattering  the 
molecule  remains  in  its  ground  electronic  level. 

Of  particular  interest  is  the  effect  on  the  cross  section  of  the  molecular  rota- 
tion. It  is  instructive  to  view  the  problem  classically,  although  the  calculation  is 
quantum  mechanical.  If  there  were  no  rotation,  the  scattering  cross  section  would  be 
given  (in  appropriate  units)  by  the  magnitude  squared  of 

f(K)  = /elK  r p(r  ) dr 

where  the  integration  is  over  all  space;  K >s  the  scattering  vector,  the  difference  be- 
tween the  wave  vectors  of  the  incident  and  scattered  neutrons.  p(r),  the  density  of 
electrons  of  unpaired  spin^  of  the  molecule,  depends  on  the  orientation  of  the  mole- 
cule. We  consider  now  two  limiting  cases  of  the  scattering  from  a rotating  molecule 
depending  on  the  ratio  a of  the  time  of  transit  of  the  neutron  across  the  molecule  to 
the  period  of  molecular  rotation.  For  a « 1^  each  neutron  sees  an  instantaneous 
picture  of  the  molecule,  so  that  the  cross  section  is  given  by  averaging  |f|  over  all  orienta- 
tions of  the  molecule.  This  limiting  case  corresponds  to  neutrons  of  high  energy.  For 
o » 1,  corresponding  to  very  slow  neutrons,  each  neutron  sees  the  same  smeared  out 
charge  distribution  given  by  averaging  p over  all  orientations,  and  the  cross  section 
is  given  by  squaring  the  corresponding  f.  In  brief,  the  cross  sections  are  given  in  the 
two  limiting  cases  a « 1 and  a » 1 by  <jf|Z>  and|<f>|2,  respectively.  The  general 
quantum  mechanical  expression  for  the  cross  section  reduces  to  the  above  expressions 
for  molecules  initially  in  J = 0 rotational  levels. 

Under  experimental  conditions  a does  not  approach  closely  either  limiting 
case,  so  that  an  investigation  of  the  scattering  for  the  intermediate  range  of  a is  in 
order.  Ari  estimate  of  the  rotational  period  may  be  obtained  by  equating  the  mean  ro- 
tational energy  of  to  kT.  For  molecules  at  room  temperature  and  incident  neutrons 
of  wave  length  5 A.  , a is  about  l/4. 

A general  quantum  mechanical  expression  for  the  cross  section  which  takes 
account  of  the  rotational  effects  has  been  obtained  in  terms  of  a number  of  integrals 
over  angle  involving  f.  The  evaluation  of  these  integrals  is  continuing. 
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15.  TRANSITION  PROBABILITIES  FOR  X-RAY  EMISSION  IN  THE  TRANSITION  METALS 


In  the  previous  Progress  Report  we  set  up  the  4s  wave  function  as  an  or- 
thogonalized  plane  wave,  which  is  by  definition  a plane  wave  made  orthogonal  to  the 
ion-core  Bloch  wave  functions  by  the  Schmidt  orthogonalization  method.  The  OPW 
corresponding  to  the  4s  crystal  wave  function  is 


x(r ) 


Z ^ +j(ir) 


h = 


/• 


ik. 


<}>.*(r ) dT 

J - 


J 4>j*(r  ) 4>j(r  ) dT 


The  coefficients  involve  the  Fourier  transform  of  the  Hartree  atomic  wave  functions 
4»j(r ).  These  are  now  being  calculated  by  expanding  the  plane  wave  in  terms  of  spheri- 
cal waves  and  carrying  through  the  integrations. 


A.  J.  Freeman 
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16  DECOMPOSITION  OF  DIRECT  PRODUCTS  OF  SPHERICAL  HARMONICS 


Most  problems  involving  one  or  more  centers  of  force,  as  for  example  mo- 
lecular or  crystal  problems,  are  simplified  by  the  introduction  of  spherical  harmonics. 
These  functions  form  a complete  set  for  the  surface  of  a sphere,  and  so  any  function 
of  the  type  likely  to  be  met  in  physics  can  be  expanded  in  terms  of  them.  In  particular, 
the  product  of  two  or  more  spherical  harmonics  is  so  expansible.  The  expansion  co- 
efficients are  certain  integrals  which  are  simple  in  principle  but  exceedingly  tedious 
to  perform.  In  practice,  operator  methods  can  be  devised  which  depend  on  the  "up- 
ward shift"  and  "downward  shift"  operators  L+  and  L , to  be  defined  below,  together 
with  the  decomposition  of  simple  products.  Strictly  speaking,  the  only  product  which 
needs  to  be  defined  is  that  of  the  general  harmonic  with  the  three  harmonics  of  degree 
unity,  for  ail  higher  products  can  be  built  on  these  as  a basis,  but  it  is  sometimes 
convenient  to  have  explicit  expressions  for  some  of  the  higher  products.  In  the  work 
below  we  give  general  expressions  for  all  products  involving  harmonics  of  degree 
zero  through  three  inclusive,  and  specific  form  to  all  products  of  combined  degree 
not  exceeding  six. 

It  will  be  convenient  to  use  functions  of  norm  4ir,  i.  e.  , 

/e.  eT]  , d£2  = 4tt  6fJ1  6 , 

#m  I'm'  It'  mm' 

The  relationship  between  partners  of  the  same  degree  is  given  with  the  Dirac  phase 
convention 


L e,  = (L  + iL  ) ef 
im  x y'  im 


L e,  = (L  - iL  ) e, 
Im  x y'  Im 


= + 1 Ctn  8 eim  = >/(l-m)(i+m  + 1)  e,  m + l 

= e_1<1^  + ictn0i4}  e!m=y(i+m)(i-m  + ei,  m-1 


These  lead  to  the  general  relationship  e^  m = (-  l)m  e^  which  can  be  used  to  elimi- 

nate about  half  of  all  possible  products.  It  is  easy  to  show  that  an  explicit  form  for 
these  functions  is 


m 


_ /(2!  + l)(l  - m)J 

■ V it  x 


..  /U*  + 1) 

**  (I  + 


)(!  - m)! 


(I  + m).1 


(p)  ‘ (-sin  6 e^)m 
P|(m^  (p)  ‘ (sin  6 e ^J111 


m ^ 0 


m ^ 0 


/ \ ^ tl"l 

where  P^  (p)  is  the  m derivative  of  the  I Legendre  function  Pj(p)  with  respect 
to  its  argument  p = cos  6.  As  special  cases  we  have 
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■Jll  + 1 Pj(cos  8) 

^zni(.sin9e‘v 

Zl  1! 

The  expressions  which  follow  were  derived  in  at  least  two  of  the  following 

ways: 

(1)  successive  synthetic  division  of 

min(l,  i1) 

elm  el'm'  ” Y cp  ei  + f'  - 2p,  m + m' 

p = 0 

by  e ..  „ . . to  determine  the  coefficients  c . 

ill'-  2p,  m + m p 

(2)  expansion  of  e „ e.,  as  above,  followed  by  successive  application  of  L to 

XX  X f “ x ^ 

produce  e^  e^,  ^ and  then  of  L to  produce  simple  algebraic  equations  easily 

solvable  cnainwise  for  efI>  _ Jt  + j,  eJm  eJS  _ Js  + 2,  ... 

(3)  repeated  multiplication  by  the  harmonics  of  low  degree,  as  illustrated  by  the  fol- 
lowing example  (for  degree  unity):  from  e.,gl  = (e.  . / we  get  upon  operation  with 

_ ||  _ ^ I X X x 1 I 

(L  ) a homogeneous  polynomial  of  degree  i1  in  e^,  e1Q,  ej  _ ^ for  e^,m,  The 

product  ejm  can  then  be  formed  by  successive  application  of  the  associative  law 

of  multiplication  in  such  a way  that  one  factor  is  always  a harmonic  of  degree  one. 

The  expressions  were  checked  with  the  aiu  of  certain  summation  formulas 

such  as 


e00  = 1 
eio  = 

eif  = 


Y 


m = - 


i 


efm  e!m 


21+1. 


General  Formulas 


eim  e00  elm 


/ 3 {I-mfZJJ  / 3 CTTmJT  o 

elmel  - 1 V 2(21+  l)(2i+  3)  (i  - m)!  ef+l,m-l  V2(21  - 1)(2I+  1)  (l+  m - 2)!  1‘1>  m_ 

/ 3 (1  + m + l)(l  - m + 1)  r j.  / 3 ji+m)(i-m)  e 

elme10  ~v  (21+  i)(2i+  3)  1 i+l,m  V (2J  - l)(2i+ 1)  1 i_1 
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1 / 3 5 (x  - m + 4).' 

eime2  - 2 " 2(ii  + 3)  V 2(2I  + i)(2l  + 5)  (j  . m).  eJ+2fm-2 

1 /3  • 5 (I  + m)'  (I  - m + 2).' 

(21  - 1 ) (21  + 3) v 2(1  + m - 2)V  (I  - m).'  e*.  m - 2 


, 1 / 3-5  (i+m).' 

2(21  - l)V  2(21  - 3)(2i  + I)  (i+m  - 4):  6i  - 2,  m - 2 


eime2  - 1 


1 / 3 • 5 (i  + m + 1)  (i  - m + 3).' 

“ 2 i + 3 » 2(2i  + l)(2i  + $)  1 (I  - m)! 


"1  + 2,  m - 1 


2m  - 


(2i  - l)(2i  + 3) 


11  e(».) 


35 


(i  + m)!  (1  - m) 


(21  - 1)V  2(21  - 3)(2i  + 1)  (i  + m - 3).'  1 i-2,  m-1 


eime20 


- 3 / 

2(2i  + 3)V  (21  + 


(i  + m + 2).1  (i  - m + 2)? 


l)(2l+5)  (I+m).'  (I  - m).'  i+2,m 


+ Ml  i)  ~ 3mZ  Vi  e 

(21  - l)(2i  + 3)  V 


+ 3 / 5 (I  + m).'  (1  - m).'  ~ 

2(21  - 1) V (21  - 3)(2i+  I)  (1+  m - 2)!  (i  - m - 2).'  1"2.  m 


1 / W (i  - m + 6).'  _ 

eime3-3  4(21  + 3)(21  + 5)^  (21  + l)(2!  + 7)  (i  - m).'  ei+  3,  m - 3 

3 / W (1  + m)!  (I  - m + 4); 

4(21-  l)(2i  + 5)  V (21  + !)(?*+  3)71  + m - ?.)'  (l  - m).'  6i+  1 

3 / F7?  (1  + m);  (I  - m + 2).' 

4(21  - 3)(2i  + 3)V  (z*  - lj(2l  + 1)  (i+m  - 4).'  (i  - m).'  L~  1 

1 / 5 7 (1+m).' 

4(21  - 3)(2i  - 1)V  (21  - 5)(2i  + I)  (i  + m - 6).'  ' 3.  m " 3 


, m - 3 
, m - 3 


-53- 


(DECOMPOSITION  OF  DIRECT  PRODUCTS  OF  SPHERICAL  HARMONICS) 


e . e . 


e.  e. 


eime30 


1 

/ 2 • 3 • 5 • 

7 (1  + m + 1)  (1  - m + 5).' 

4(21  + 3)(2! 

+ ■5)^  (21  + 1)(21  + 

7)  1 (1  - 

m).' 

1+3,  m - 2 

i + 3m  - 2 

: / 2 • 3 • 5 • 

7 (1  + m)  (1  - m + 

W e 

4 (2i  - l)(2i 

+ 5)'V  (2|  + ^(21  + 

3)  1 (1  - m).' 

el  + 

1,  m - 2 

J - 3m  + 3 

i / 2 • 3 • 5 • 

7 (1  + m)I  (1  - m + 1) 

4(21  - 3)(2I  - 

1 3)  V (21  . !)(2|  + 

1)  (1  + m - 3).' 

1 

el-  1,  m -2 

1 

/2  • 3 • 5 • 

7 (1  + m)!  (l  - m) 

4(2/  - 3)(2! 

- 1)  V (21  - 5)(21  + 

1)  (1  + m - 5).'  1 

i - : 

>,  m - 2 

5 

/ 3-7 

(1  + m + 2).'  (1  - 

m + 4), 

1 

4(21  + 3)(2 1 

+ 5)  v (21  + 1)(21  + 

7)  (1  + m)!  (1 

- m)J 

cl+  3,  m - 

i2  - lOim  - 

15m2 +71+  5m  / (l  - : 

m + 2).' 

4(21  - 

1)(21+  5)  v (2|+  i)(21+  3)  (1 

- m): 

1+  1,  m - 1 

I2  + 101m  - 

15m2  -51+1  5m  - 

6 / 3-7 

(1  + m).'  „ 

4(21  - 3)(21  + 3) 

v (21  - l)(21  + 1)  (1  + m 

„ 2).'  J 

5 

/ 3 7 

(1  + m)!  (l 

- m).' 

4(21  - 3)(2! 

- l)  V(2 1 - 5)(21  + 

1)  (1  + m - 4).'  (1  - 

m - 2) 

; C1  - 3,  m - 

5 

/ 7 

(1  + m + 3).'  (1  - 

m + 3).' 

2(21  + 3)(21 

+ 5)  v (21  + 1)(21  + 

7)  (l  + m)J  (1 

- m).1 

1 + 3,  m 

3[Kl+  2)  - 

5m  “1  / 7 

(1  + m + 1 )(l  - 

m + 1) 

2(21  - 1)(21 

! + 5)  V(21  + 1)(21 

+ 3)  1 

1 

el+  1,  m 

3[(!+  1 )(1  - 

1)  ~ 5m2]  f 

7 (1  + m)  (i 

' - m) 

2(21  - 3)(21  +3)  v (21  - 

l}(21  + 1)  1 

1 C 

1 - 1,  m 

5 

/ 7 

(l  + m)!  (1 

- m).' 

2(21  - 3)(2! 

^17  V (H  - 5)(21  + 

1)  (l  + in  - 3).!  (1  - 

m - 3).' 

el-  3,  m 
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I + I'  = 2 


Specific  Formulas 


i + i'  = 3 


i + V =4 


e,,  = Vi  3/5 


Ti  ll 


6 1 1 e 1 0 


'22 


= Syl 


e 1 1 e i - 1 = ^7* 


'21 


e20  ” e00 


e 1 oe  1 0 


2 Tl/i  e,„  + e. 


'20  00 


e22el 1 


e22e 1 0 


zVTfi 


33 


= Ti/? 


'32 


e22Cl  - i = ^3/5  ‘ 7 e30  " ' 3/5  e 


10 


e 2 1 e 1 1 ~ 7 


'32 


e 2 1 e 1 0 


2-x/2  • 3/5  • 7 e31  + TV* 


li 


e2 le 1 - 1 


3 Vl/5  • 7 e30  - 7^5 


10 


e20e 1 1 


e20e  1 0 


3^2/5  • 7 e31  - y/Tfl  e 


11 


3^3/5  • 7 e30  + 27175 


10 


e 33e  1 1 


e33e  1 0 


27T73 


'44 


= 7T/3  e 


43 


e33e  1 . ! = 7l/3  • 7 e42  - 3 71/7 


22 
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I + V = 4 (con'd) 


e32e 1 1 
e32e 10 
e32el - 1 

e 3 1 e 1 1 
e 3 1 e 1 0 
6 3 1 e 1 - 1 

e30el  1 
e30e 1 0 

e22e22 
e22C21 
e22C20 
e22e2 - 1 
e22e2 - 2 


1 • c43 

2 e42 

+ v'sT?  e22 

VT/7  e41 

- 

3/ 7 e2, 

7iT7e42 

- 7v ? 

• 7 e 
e22 

'/ST?  e4l 

+ 2 7T"- 

3/5-  7 ®21 

e4Q 

- 3^/2/ 

5 • 7 "20 

72  • 5/3 

?e41  ' 

' 3^175"  7ez, 

4 7l/3 

7 e40  + 

373/5  ? e?0 

v^e43 

(1/7)  yr~5  e42  - (2/7)7?  e22 
(i/7)y5e41  - (1/7)  ^2  • 3 • 5 e21 

( 1 / 7)  e40  " e20  + 1 e00 


e21e21  = (2/7)72  • 5 e42  + (l/7)  72-3-5  e22 
e21e20  = (1/7)  72  3 • 5 e41  + (1/7)7?  e21 
e21e2  - 1 = ^4/7)  e40  ” Cl/7)75’  e2Q  - 1 eQ0 


'20  20 


= *6/7)  s„n  + (2/71  ./?  p.  . + 
7 7 40  7 v " 


1 

uu 


e44el  1 = v^3  • 5/11  e55 

e44e10  = ^3/  1 1 e54 

e44er.  1 = v'1/3  • 11  e53  - 27T/1  e33 
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I + I1  = 5 (con'd) 

e43e i j — 2 73/11  654 

e43ei0  = 4 7l  / 3 i~i  e53  + 7l/  3 e33 
e43e  1 - 1 = TV11  e52  " 1 e32 

e42en  = 2^7/3  • lle^  - ^/l/3  • 7 e33 

e42e  10  — ”77 1 1 e52  2 ^ ®32 

e42ei  _ j “ 72/  1 1 e^j  - ^/sjl  e^j 

e41el  1 = 77/11  e52  ‘ 71/7  e32 

e41e10  = ^ 72/ 1 ^ e5j  V 5/7  e^j 

e41el  . i = V2  • 5/3  • 11  e5Q  - 72  5/3-7  e30 

e40el  1 = v/^_1  e5 1 - T2/7  e31 

e40e10  = 57l/3  • 11  e5Q  + 47l/3  • 7 e,0 

e33e22  = 572/3  1 1 e55 

e33e21  = 2 75/ 3 11  e&4 

e33e20  = (2/3)  75/11  e53  - (l/3)75e33 

e33e2  - 1 = (l/3)  Ts/Tl  e52  ‘ <73)75  e32 

e33e2  - 2 = (1/3)75/7  • 11  e51  - (1/3)7^  e31  + 

e32e22  = 72 ' 5/11  e54 

e32e2 1 = <4/3)  75/11  e53  + (73)75  e33 

e32e20  = ^5/ 1 1 e52  + 0 

e32e2-  1 = 272  - 5/3  • 7 11  e51  - 773  e31  - 72  3/7  en 
e32e2-2  = (5/3)71/7-  11  e60  - (2/3)  e3Q  + 73/7  e10 
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I + I1  = 5 (con'd) 

e31e22  = (5/3)  e53  - (1/3)72  e33 

e31e21  = (5/3)-y37Tl  e52  + Jl/3  e32 

e31e20  = 5 ' 11  e5i  + e31  + 372/5  • 7en 

e 3 1 e 2 _ 1 = (5/3)72 ' 5/7  • 7l  e5()  - e3Q  - 2^2  -3/5-7  e1Q 

e 3 1 e2  _ 2 = 5 7l/  3 • 7 • 11  e5  _ j - Zjl/l  • 5 e3  _ j + 73/5  • 7 ej  _ j 

e30e22  - (^/^)  e52  ” (^/3)  e32 

e30e2 1 = (20/3)yi/7  • 11  e51  + (1/3)75/5  e3J  - 37l/5  • 7 en 

e30e20  = (10/3)757?  • lle50  + (4/3)7I75e30  + 373/5  • 7 e1Q 


55el  1 

= 

3^2/13  e66 

55C 1 0 

= 

'/3/13e65 

55e  1 - 

1 = 

73/11  • 13  € 

!64  ‘ 

■ 7T 

• 5/11  e44 

54ell 

= 

73  • 5/13  e65 

54e  1 0 

— 

273  • 5/11  • 

13 

e64  + 

"'^3/ 1 i e44 

54el  - 

I = 

37l/ll  • 13 

e63 

- 2 73/n  e43 

53e  1 1 

— 

373  • 5/11  • 

13 

e64  " 

7l  / 3 • lie 

5 3e  1 0 

= 

97l/ll  13 

e63 

+ 471/3  • 11  e43 

53e  1 - 

1 = 

372/11  • 13 

e62 

- 27 7/3  • 11  e42 

526 1 1 

= 

6 73/11  • 13 

e63 

- Tl/n  e43 

52e  1 0 

= 

472  • 3/11  • 

13  1 

e62  + 

^7/11  e4i 
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e52el  - 1 


'Z  3 • 5/11-  13  e,  . - ./  7 / 1 1 


e51eu  = 2^3  • 7/ 1 1 13  eA?  - yi/Tl 


e 5 1 e 1 0 
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e50el 1 
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“43e22 

3 P 

“43e2 1 
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e42e2-2  = 5/11  yi/TT  e6Q  - 2 • 3 • 5/7  • ll-TIe^  + l/?y3  • 5 e20 


(DECOMPOSITION  OF  DIRECT  PRODUCTS  OF  SPHERICAL  HARMONICS) 

I + V = 6 (con'd) 

e41e22  = 3/11-/2  • 5 • 7/13  e63  - 3/11^2  -3-5/7  e 43 

e4le21  = 4/11^2  • 5 • 7/13  e&2  + 9/7  • - 1/7/?  e22 

e41e20  = 5/llVz  ' 3 ' 7/13  e6l  + 17/7  ' H-/5e41  + l/7^—3_5e21 

e41e2  _ ! = 2 • 5/1172-  3/T3  e60  - 5/7  • 11^2  • 3 e4Q  - 1/7^2  -3-5  e2Q 

e41e2-2  = 5/  1 i 77713  e6  _ 1 - 2 5/7-  1172  -_3  • 5 e4  _ j + 1/7^5  e,  , 

e40e22  = 72  7 / 1 5 e^2  “ 2 - 3 - 5/  < • 11 7^  e4£  ^ / 7 e22 

e40e21  = 5/1175  • 7/13  e6l  + 5/7-  ll72  7 3 e4l  - 4/7e2[ 

e40e20  * 3 • 5/"'/s/T3  e6Q  t 20/7  • 11 7?  t 2 • 3/7  e2Q 

e33e33  = 2 ‘ 577/3  • 1 1 ” ’5e66 

e33e32  = 572  • 7/3  • 11*  13 

e33e31  = 5/1172  • 5 • 7/3  • 13  e64  - 1/ 1 1 72  - 3 • 7 e44 

e33e30  = 2 " 3/ll77/3  l"3  - 3/ 1 1 T7"  e43 

e33e3  - 1 = 2 ‘ 5/3  ‘ 11  -n/7/^  e62  - 3/1172  • 3 e42  + 1/372  e22 

e33e3  - 2 = 5/3  ’ 1 17 7/l3  e6l  - 1/  1 1 e41  + 1/375  e2i 

e33e3-3=5/3  11  ^l/  I3  e60  ” 3/  1 1 e40  + l/37^  e2Q  - 1 eQ0 

e32e32  = 2 * 5/1172  • 7/13  e64  + l/U^/2  -5*7  e44 

e32e31  = 5/1172  • 3 • 7/13  e63  + l/ll72  • 7 e43 

e32e30  = 4‘  5/3  ‘ 1172 ' 7/l3  efe2  - i/ 1 1 T3  e4^  - 2/3e22 

e„e,_  j = 5/1175  • 7/3  • 13  e6l  - 4/llTF e4J  - 71/1  e21 
e32e3  -2  = 2 5/11  7l / 1 3 e6Q  - 7/ll  e4Q  + 0 e2Q  + 1 eQ() 


-60- 


(DECOMPOSITION  OF  DIRECT  PRODUCTS  OF  SPHERICAL  HARMONICS) 


I + V = 6 (con'd) 

e31e31  = 2 * 5/ll75  • 7/3  • 13  + 2/11^2  • 5 e42  + 2^2/3  • 5 e22 

e31e30  = 5 ' 5/3  • 11^2  • 7/13  e6]  * e41  + l/lVZJl  e£1 

e 3 1 e 3 _ i = 5 5/ 1 1 VT/Ti  e^0  - 1/  1 1 e4Q  - a/T/5  e2Q  - 1 eQ0 

e30e30  = 4 * 5 5/  3 • i 1 -</l  / 1 3 e^Q  + 2 • 3/ 1 1 e4Q  + e20  + 1 eQ0 

H.  C.  Schweinler 
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17,  THERMAL  VIBRATIONS  OF  Cu-Zn  SYSTEM  CRYSTALS 


Work  on  this  topic  continues  to  be  focused  almost  entirely  on  the  determina- 
tion of  the  force  constants  between  the  various  pairs  of  atoms  in  each  crystal.  ^ 
Workers  in  Professor  B.  E.  Warren's  x-ray  laboratory  at  M.  I.  T,  have  succeeded  in 
working  backward  from  their  thermal-diffuse -scattering  isodiffusion  contours  in  re- 
ciprocal lattice  space  to  obtain  values  for  the  interatomic  constants  by  finding  the 
best  match  possible  between  the  theoretical  and  experimental  intensity  curves  along 
certain  simple  directions.  ^ This  technique  furnishes  the  only  experimental  values 
available  at  present  for  atomic  force  constants.  However,  the  experimental  measure- 
ments are  quite  difficult  because  of  the  low  intensities  and  possible  surface  effects 
involved.  Likewise,  higher-order  scattering  terms  and  higher  than,  say,  third-nearest- 
neighbor  interactions  cannot  be  accurately  taken  into  account  by  the  theory  used  be- 
cause of  the  involved  nature  of  the  calculations;  nor  can  the  Compton  background 
scattering  be  corrected  accurately  everywhere.  It  is  important  to  have  theoretical 
atomic  constant  values  for  comparison  with  the  x-ray  values  because  of  these  weak- 
nesses as  well  as  for  their  own  intrinsic  interest. 

Only  elastic  constant  calculations  have  so  far  been  published,  and  they  are 
quite  different  from  atomic  constant  calculations  in  that  elastic  deformations  are 
homogeneous,  while  the  atomic  displacements  to  be  considered  are  localized  ones, 
much  as  in  typical  impurity  problems.  The  procedure  (analogous  to  that  in  the  el- 
astic case)  of  taking  the  second  mixed  partial  derivatives  with  respect  to  two  nuclear 
coordinates  of  the  crystal  cohesive  energy  to  obtain  the  atomic  force  constants  would 
consequently  be  a difficult  task.  Fortunately,  Feynman's  theorem  tells  us  that  the 
force  on  one  nucleus  at  its  lattice  site  resulting  from  the  displacement  of  another  in 
an  otherwise  perfect  lattice  is  the  classical  Coulomb  force  produced  by  the  change  in 
total  charge  distribution  of  all  the  electrons  in  the  crystal,  bound  and  free,  from  the 
equilibrium  distribution  which  produces  no  net  forces.  ^ Since  the  limiting  value  of 
such  a force  as  the  displacement  goes  to  zero  is  in  fact  just  the  atomic  force  constant, 
the  calculation  of  the  perturbation  in  electronic  wave  functions  resulting  from  a small 
nuclear  displacement  is  the  only  difficult  step  in  computing  an  atomic  constant  by  this 
Feynman  approach,  which  was  suggested  by  Professor  Slater.  Moreover,  the  force 
constants  for  all  orders  of  neighbors  can  be  easily  obtained  by  classical  Coulomb  po- 
tential integrations  from  this  one  expression  for  the  change  in  charge  density. 

The  difference  between  the  perturbed  and  original  nuclear  frameworks  of 
charge,  which  is  the  perturbation  acting  on  the  electron  cloud,  is  clearly  a dipole  at 
the  site  of  the  displaced  nucleus.  Three  lines  of  attack  have  been  followed  out.  Treat- 
ing the  conduction  electrons  as  a Thomas -Fermi  gas  spread  out  through  the  whole 
crystal,  one  obtains  by  a standard  approximation  the  result  that  the  ion-core  dipole 
as  screened  by  the  electron  cloud  has  a charge  distribution  and  a potential  which  are 
both  proportional  to 
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(pcos_e)(l  + qr)  qr 

This  potential  gives  considerably  more  anisotropic  force  constant  values  than  a pure 
dipole.  Furthermore,  for  a typical  metal  the  parameter  q is  such  that  the  screening 
takes  place  within  a short  distance,  and  in  fact  is  about  90%  complete  within  a radius 
equal  to  that  of  the  actual  ionic  core.  It  would  follow  that  the  screened  dipole  force 
at  a first -neighbor  nucleus  would  only  be  a few  percent  of  that  produced  by  the  un- 
shielded ion-core  displacement.  Clearly  these  results  cannot  be  relied  upon  quan- 
titatively, since  the  conduction  electrons  are  actually  not  free  to  move  in  the  ion-core 
volumes  in  which  the  model  of  Thomas -Fermi  has  them  doing  most  of  the  screening. 
Next  an  attempt  was  made  to  take  the  second  step  in  a proper  self-consistent  field 
procedure  by  using  the  Thomas-Fermi  screened-dipole  potential  as  the  perturbation 
term  in  a quantum -mechanical  treatment.  However,  even  using  the  simplest  free- 
electron  model  of  a metal  and  the  equivalent  of  the  Born  scattering  approximation  in 
solving  the  one -electron  Schrodinger  equation,  the  integral  for  the  perturbed  wave 
function  was  quite  difficult.  Since  further  integrations  over  all  occupied  free -electron 
states  and  then  over  all  space  with  the  Coulomb  factor  l/r  would  have  been  required 
to  get  the  new  screened  potential,  it  was  not  worthwhile  to  carry  out  the  calculation 
in  view  of  the  crudeness  even  of  a free  electron  model  which  excludes  the  conduction 
electrons  from  the  central  ion  core. 

A number  of  general  statements  cam  be  made  regarding  the  atomic  force  con- 
stants. From  the  equations  of  motion  of  the  thermal  vibrations  and  the  symmetry 
properties  of  the  lattice  one  cam  easily  derive  symmetry  relations  for  the  force  con- 
stants, and  this  has  been  done  for  the  first  three  orders  of  neighbors  in  the  fee  and 
bcc  lattices.  Three-by-three  matrices  contain  all  the  force  constants  between  a 
given  pair  of  atoms,  and  the  symmetry  restricts  the  internal  structure  of  these  ma- 
trices as  well  as  relating  the  matrices  for  equivalent  pairs  of  atoms.  For  example, 
there  turn  out  to  be  but  four  independent  parameters  in  the  twenty-four  three -by -three 
force  constant  matrices  relating  the  (arbitrary)  origin  atom  to  the  twenty -four  third 
nearest  neighbors  in  the  fee  lattice.  Elastic  constants  for  the  CuZn  system  crystals 
have  been  quite  accurately  measured  experimentally,  although  the  theoretical  values 
are  not  very  reliable.  Thermal  vibrations  include  the  long  wavelength  elastic  vibra- 
tions as  a special  case,  and,  as  one  would  expect,  the  elastic  constants  can  be  ex- 
pressed as  sums  of  combinations  of  the  atomic  constants  for  each  order  of  neighbors. 
Born  has  proved  this  general  result  by  comparing  the  leading  terms  in  a perturbation- 
theory  expansion  of  the  normal  coordinate  equations  of  motion  for  thermal  vibrations 

(with  the  wave  vector  as  perturbation  parameter)  with  the  terms  of  the  elastic  wave 

(4) 

equation.  ' These  relations  have  been  established  explicitly  through  third  neighbors 
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in  the  bcc  and  fee  lattices  in  the  course  of  our  work.  They  are  useful  as  a completely 
unambiguous  partial  check  of  the  atomic  constants  with  reliable  experimental  results, 
although  of  course  there  are  many  more  independent  force  constants  than  elastic  con- 
stants even  in  the  first  few  rings  of  neighbor  atoms.  Born  lias  further  proved  by  a 
thermodynamic  argument  that  a lattice  satisfying  the  usual  elastic  or  long  wavelength 
criteria  for  stability  also  automatically  is  stable  to  all  the  thermal  vibrations:  i.  e. , 

7 7 -a 

that  w > 0 for  all  the  10  or  so  normal  vibrations  if  only  the  few  elastic  stability 
criteria  are  satisfied.  ^ Thus  agreement  of  the  appropriate  sums  of  atomic  constants 
with  the  experimental  elastic  constants  guarantees  that  all  stability  criteria  are  satis- 
fied as  well,  although  there  are  still  a few  microscopic  criteria  which  it  proves  in- 
structive to  evaluate.  It  should  further  be  mentioned  that  the  order  of  magnitude  of 
the  observed  elastic  constants  shows  that  the  actual  atomic  constants  are  of  the  same 
general  magnitude  as  those  which  would  result  from  an  unscreened  ion-core  displace- 
ment, contrary  to  the  Thomas-Fermi  model  prediction. 

At  present  a general  semi-empirical  study  of  the  force  constants  is  being 
carried  out.  From  the  Thomas-Fermi  model  results  one  feels  fairly  confident  that 
most  of  the  perturbation  in  electronic  charge  distribution  takes  place  within  the  sphere 
through  the  nearest  neighbor  atoms.  Hence  this  Ap  can  be  approximated  in  the  usual 
way  by  a multipole  expansion  valid  at  all  neighboring  atoms.  The  force  constant  ma- 
trices produced  by  each  multipole  for  the  first  through  third  orders  of  neighbors  have 
been  computed  for  fee  and  bcc  lattices  in  terms  of  the  multiple  moments  as  parameters. 
These  multipole  moments  are  then  chosen  so  that  the  sum  of  the  contributions  they 
make  to  the  force  constants  relating  various  orders  of  neighbors  agree  with  the  values 
which  give  the  best  fit  with  the  x-ray  data.  Finally  the  Ap  can  be  built  up  from  these 
multipole  moments.  Since  the  multipole  potentials  each  have  a different  radial  de- 
pendence, the  force  constants  produced  by  each  multipole  at  all  orders  of  neighbors 
must  separately  satisfy  the  symmetry  requirements.  In  the  fee  and  bcc  structures 
symmetry  thus  rules  out  monopole,  quadrupole,  etc.  terms.  Various  other  more 
qualitative  restrictions  on  and  relations  among  the  force  constants  have  been  found 
and  can  be  applied  in  determining  a satisfactory  match.  In  bcc  a-Fe  fairly  good  match- 
ing with  x-ray ^ and  elastic  values  has  been  obtained  for  first  and  second  neighbors  using 
dipole  plus  octopole  terms.  The  x -ray  atomic  force  constant  values  for  fee  Cu  and  hp 
Zn  being  obtained  by  Warren's  group  are  not  yet  well  enough  established  to  permit  a 
reliable  multipole  fitting;  however,  from  the  stability  criteria  it  is  clear  that  dipole 
and  octopole  terms  alone  cannot  be  enough.  At  distant  neighbors  only  dipole  terms 
will  be  appreciable,  and  it  therefore  may  prove  feasible  to  include  the  effect  on  the 
vibration  frequency  spectrum  of  all  the  further  neighbor  lattice  interactions,  using 
the  techniques  applied  in  ionic  lattices  and  the  dipole  moment  which  best  fits  the  x-ray 
constants  for  the  first  few  rings  of  neighbors.  However,  it  is  anticipated  that  the 
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primary  application  of  these  semi -empirical  results  will  be  in  forming  a clear  physi- 
cal picture  of  the  charge  distortion  resulting  from  a nuclear  displacement  as  a guide 
in  the  fundamental  task  of  developing  an  accurate  quantum -mechanical  calculation  of 
this  charge  distortion,  and  thence  the  force  constants,  from  first  principles:  several 
possible  approaches  to  this  task  are  under  consideration. 
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